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Kahler metric, with oblique (non-commuting) fluxes in magnetized brane construc- 
tions is an interesting unresolved issue, in view of the importance of such fluxes for 
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with the complex structure of the torus. By evaluating the overlap integral of the 
wavefunctions, wc give the expressions for Yukawa couplings among chiral multiplets 
arising from an arbitrary set of branes (or their orientifold images). The method 
is based on constructing certain mathematical identities for general Riemann theta 
functions with matrix valued modular parameter. We briefly discuss an application 
of the result, for the mass generation of non-chiral fermions, in the SU{5) GUT 
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1 Introduction 



Close form expressions for Yukawa couplings have been written down for string construc- 
tions involving branes at angles [1,2] or those with magnetized branes [3-17]. In the IIA 
picture, the interaction is described by the worldsheet instanton contributions from the 
sum of areas of various triangles that are formed by three D6 branes intersecting at three 
vertices, forming a triangle. This is due to the fact that the intersection of branes relevant 
for Yukawa interactions are those which are point-like giving chiral multiplets. Line or 
surface like intersections, on the other hand, would give rise to interactions of non-chiral 
matter. In these discussions, the orientation of the branes themselves are parameterized 
by three angles in the three orthogonal 2-planes, inside T^. These results have been fur- 
ther generalized to include Euclidean D2 brane instanton contributions to the Yukawa 
couplings [18-26], generating up quark and right handed neutrino masses through a Higgs 
mechanism, in a particular class of models. A limitation on the exercise performed in these 
papers comes from the factorized structure of the tori, which arises from the orientations of 
the brane wrappings that are classified by angles in three different planes, rather than 
their general orientations in the internal six dimensional space parameterized for instance 
by the SU{3) angles in super symmetric situations. 

Similar results for perturbative Yukawa couplings have also been obtained in the mag- 
netized brane picture, based on their gauge theoretic representation [5]. In this case, the 
interactions are given by the overlap integral of three wavefunctions (contributing to the 
interaction) along internal directions. The wavefunctions correspond, in the ordinary field 
theory context, to those belonging to two fermions and a scalar, and are given by Jacobi 
theta functions, when fiuxes are turned on along three diagonal 2-tori. The relationship 
between the Yukawa interactions in the magnetized brane constructions and those involv- 
ing D6 branes, have also been estabhshed using T-duality rules. However, these exercises 
have once again been of limited scope due to the fact that explicit expressions are written 
down only for magnetized branes with fluxes that are diagonal along three T^'s. 

Technically, the wavefunctions of chiral fields participating in Yukawa interactions are 
defined in terms of Jacobi theta functions, with a modular parameter identified as a 
product of the complex structure of the T^, with the fiux that is turned on along it. The 
Yukawa interactions are therefore computed for the case when the six dimensional internal 
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space is of a factorized form: 

T'^ xT'^ xT'^ e T^. (1.1) 

For applications to model building with moduli stabilization though, one in general 
needs to include both 'diagonal' and 'oblique' fluxes [10,12,13,27]. Therefore it is imper- 
ative that we generalize previous results further and obtain interactions involving branes 
with oblique fluxes. As stated, in the language of D6 branes such generalizations would 
amount to intersections of branes with orientations given by SU (3) rotation angles, re- 
sulting to = 1 supersymmetry in D = 4 with chiral matter. In view of the importance 
of such fluxes in obtaining realistic particle physics models with stabilized moduli, and to 
describe the interactions among the chiral flelds, we shall study the explicit construction 
of fermion (and scalar) wavef unctions on compact toroidal spaces with arbitrary constant 
fluxes. 

Scattered results on fermion wavefunctions in presence of constant gauge fluxes, on tori 
of arbitrary dimensions, exist already in the literature [5,28]. However, they are of limited 
use for our purpose. First, any wavefunction obtained through a diagonalization process 
of the gauge fluxes [28], is not in general suitable for obtaining an overlap integral of 
wavefunctions. This is because the flux matrices need not commute along different stacks 
of branes that participate in the interaction through the chiral multiplets, arising from the 
strings that join these branes and therefore they are not simultaneously diagonalizable. 

In [5], a set of wavefunctions was given for constant gauge fluxes. However, once again, 
explicit results are valid only for those fluxes which satisfy a set of 'Riemann conditions', 
including a positivity criterion on the flux matrices. As the analysis in our paper will 
clarify, the positivity restrictions on the fluxes is due to the fact that the given wavefunction 
in [5] corresponds to a speciflc component of the 2" dimensional Dirac spinor for a 2n- 
dimensional torus T^". We will show that this restriction is relaxed, if one considers 
wavefunctions of various chiralities, such that all possible flux matrices are allowed, though 
in our case we restrict to only those fluxes that are consistent with the requirements of 
space-time supersymmetry . 

In fact, we give explicit solutions for the wavefunctions for arbitrary fluxes, that are well 
defined globally on the toroidal space. We also give explicit mappings among the wave- 
functions of different chiralities, satisfying different consistency criterion. These mappings 
are shown to relate wavefunctions corresponding to different fluxes and complex structures 
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of the tori. We in fact reconfirm that our wavef unctions, as well as mappings are indeed 
correct, by showing that equations of motion also map into each other for the fermion 
wavef unctions just described, corresponding to different internal chiralities. 

Apart from the lack of enough knowledge about the fermion wavef unctions, the limi- 
tations on available information about the Yukawa couplings for general gauge fluxes also 
arose from the technicalities in dealing with general Riemann theta functions that are 
used for defining the wavefunctions on toroidal spaces. Internal wavefunctions of chiral 
fermions participating in the interaction are given by a general Riemann theta function 
whose modular parameter argument is determined in terms of the complex structure of 
as well as the 'oblique' fluxes that we turn on. Hence, the limitations on available results 
for Yukawa interactions in the literature, arise due to the intricacies involved in evaluat- 
ing the overlap integrals of the trilinear product of general Riemann theta functions over 
the six dimensional internal space. In particular, even for positive chirality wavefunctions 
along the internal given in [5], one finds that theta identities [29] need to be further 
generalized, in order to compute the Yukawa interactions with oblique fluxes. The task 
goes beyond the identity given in [29], since one needs to evaluate the overlap integral of 
three wavefunctions, all having different modular parameter matrices as arguments, due 
to the presence of different fluxes along the three brane stacks involved in generating the 
Yukawa coupling. 

In this paper, first, we generalize the identities used in [5] (available from mathemat- 
ical literature [29]) for the known positive chirality wavefunctions to those with general 
Riemann theta functions representing the fermion wavefunctions. This gives an explicit 
answer for the Yukawa interaction in a close form and generalizes the results of [1,5]. In 
particular, we generalize the result further for the positive chirality wavefunction, when 
general (hermitian) fluxes with all nine parameters are turned on rather than the six 
components considered before. 

On the other hand, as already stated earlier, we give explicit construction of the other 
spinor wavefunctions, as well. In these cases too, we obtain the selection rules among 
chiral multiplets giving nonzero Yukawa couplings. Now, however, the flnal answer is 
left as a real flnite integration of a theta function, over three toroidal coordinate variables. 
The integration can be evaluated numerically for any given example. Finally, in the paper, 
we also briefly discuss the issue of mass generation for non-chiral multiplets given in the 
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SU{5) GUT model, constructed in [27]. Although the Yukawa couplings can be used for 
giving the precise masses in our model, where all close string moduli are fixed, we avoid 
entering into these details in this work. 

The plan of the rest of the paper is as follows. In Section 2, we start by discussing the 
general setup, including the gauge fiuxes that can be turned on, in a consistent manner. 
In Section 3, we review the known results on the Jacobi theta identity given in [29] and 
present a proof of its validity. We also give an expression for the Yukawa interaction for 
factorized tori and 'diagonal' fluxes using the theta identity. In Section 4, we construct 
a similar identity, now for the general Riemann theta function. We then use this new 
mathematical relation for writing down the expression for the Yukawa interaction when 
oblique fluxes are present and satisfy the 'Riemann conditions' of [5]. Results are further 
generalized to include the most general flux matrices consistent with supersymmetry and 
'Riemann condition' requirements. In order to relax the later, in Section 5, we present 
the generalizations to include the wavefunctions of the other internal chiralities, in order 
to accommodate general fluxes consistent with supersymmetry restrictions. In Section 
6, we briefly present an independent analysis of the superpotential and D-terms for the 
model of [27] in order to show how masses for several non-chiral fermion multiplets can 
be generated, without evaluating explicitly the superpotential coefficients, which would 
go beyond the scope of our present work. Conclusions are presented in Section 7. In 
Appendix A we give the chiral fermion wavefunctions in the presence of constant fiuxes. 
Appendix B contains information on fiuxes in terms of windings and Chern numbers, 
while Appendix C gives some details of our model in [27] , needed for the mass generation 
analysis of Section 6. 

2 Fluxes 

We now start by describing the gauge fiuxes that can be turned on along internal tori. A 
general gauge fiux, on with coordinates = (x*, y'^), i — 1, 2, 3, has the form: 

F = pijdX^ A dX-^ 

— Px'xidx'^ dx-' + Pyiyjdy^ A dy-' + p^iyjdx^ A dy-' + Pyi^jdy^ A dx-' . (2-1) 
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Then using the definition of a general complex structure matrix il: 

dz' ^ dx' + Q]dy^ , d? ^ dx' + n^dy^ , (2.2) 

we obtain: 

F = F,i,j dz' A dz^ + F.i^j {idz' A dz^ ) + F^i^j dz' A dz^ . (2.3) 
Choosing the basis e*-' of the cohomology H^'^ to be of the form 

e^^ ^ idz' A dz^ , (2.4) 

we obtain: 

F,i,J = (Q - Q)-l^ {n'^P^.n - Cfp,y + Pl^n + Pyy) (^ " Q)-^ (2.5) 

and 

F,i-,i = (-i)(Q - n)-^^ {n^Pxx^ - ^'^Pxy+Ply^+Pyy) - ^)~^ ■ (2.6) 

In addition, F^izi is complex conjugate to and F^i^j = —F^jgi. 

Now, supersymmetry demands all fluxes to be of (1, 1) form which gives us the condi- 
tion: 

(O^P^^O - ffp^y + Ply^ + Pyy) ^ 0, (2.7) 

or equivalently: 

(Q^P^^Q - ViFp^y + Ply^ + Pyy) = 0. (2.8) 

Eqs. (2.7) and (2.8) together give two real matrix equations. These equations can then be 
used to eliminate some of the variables and write the final (1,1) form in terms of certain 
independent variables only. 

Using eq. (2.8), eq. (2.6) reduces to the following form, 

Fz^-zi = -i {Pxx^ - Pxy) - ^)~^ . (2.9) 
On the other hand, use of eq. (2.7) in eq. (2.6) gives, 

F,,,, = - n)-'''' (-nVx ■ (2.10) 

We also notice that the (1, 1) form F^izJ given in eq. (2.6) satisfies the hermiticity property: 
Fziz^ = F^izJ- To explicitly see that, we use eqs. (2.9), (2.10). 

^ _i(Q _ fi)-i^ -pJJ = Fz.-z, (2.11) 
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There are some special cases, however, in which eqs. (2.7) and (2.8) simphfy further and 
the resulting F^i^j can be written more compactly. One such case arises when p^x and Pyy 
components are turned off. In such a situation F(2,o) — condition (2.8), reduces to: 



(2.12) 



Thus far, we have concentrated on the spatial components of the gauge fluxes, but 
ignored the gauge indices. In the magnetized D-brane construction, gauge quantum num- 
bers arise from the Chan-Paton factors associated with the end points of the open strings 
for a given stack of branes. The simplest possibility is to consider fluxes with gauge indices 
given by an n X n identity matrix for a stack of £)-branes: 



(2.13) 



with m an arbitrary integer giving the 1st Chern number. All spatial indices of the gauge 
flux above have been suppressed, which are given as in eq. (2.1) by the components : Pxiyj, 
Pxixi-i Pyiyi-i whilc their quantized values are expressed in eq. (B.4) for general wrapping 
of the branes. Actually, eq. (2.13) corresponds to the situation when all the wrapping 
numbers are trivial: n^' = n^' = 1 in the denominators of eq. (B.4). F, then represents a 
stack of n magnetized D-branes with a C/(l)" gauge flux. The first Chern number for each 
of the U{1) fluxes is equal to m. Moreover, D-brane wrapping numbers on the internal 
directions, are all unity, given by a diagonal embedding of the brane in target space and 
winding around each 1-cycle once. In most of the paper, we will consider fluxes of the 
above type. 

For multiple stacks of rii branes with respective 1st Chern numbers m^, the flux matrix 
is of block diagonal form: 



m2/, 



n2 



(2.14) 



V 



•p Tip j 



and corresponds to gauge fluxes in the diagonal [/(l)'s of U{n\) x Uiu'i) x ■■■ gauge 
group. Chiral fermion wavef unctions, on the other hand, belong to the bifundamental 
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representations {ua, rib) etc. and are determined by the field equations in terms of the 
difference of fluxes in the respective stacks, parameterized by a matrix N: 

^^=Pty-ply (2.15) 

Using (2.12), it foUows that: 

(Nnf = (NQ). (2.16) 

Moreover, further conditions are imposed on the matrix N in order to have weU defined 
bifundamental wavef unctions. These are the so-caUed Riemann conditions [5] and are 
written as: 

Nj, e z , 

N./mO > 0. (2.17) 

The first condition in eq. (2.17) is the integrahty of the elements of N, that we discuss 
later on, in the absence of any non-abelian Wilson lines [5], following from the Dirac 
quantization of fluxes. To understand the last condition of eq. (2.17), one rewrites the 
(1, 1) form F^i^j, for the case when p^x — Pyy = 0. Indeed using eq. (2.12), one obtains: 

F,i,j^-iPxy{n-n)-\ (2.18) 

which matches with the expression for H in eq. (4.73) of [5] upon the identiflcation 
N"^ = Pxy and H — ^N'^.Imfl"^. The positivity requirement on H then arises from the 
condition that the solutions of the Dirac equation, corresponding to chiral wavef untions, 
be normalizable. 

Gauge fluxes on branes with higher wrapping numbers can also be given a gauge 
theoretic interpretation. The method, as stated earlier, is based on a representation of 
the magnetized brane constructions [5] in terms of fluxes along internal directions in a 
compactifled gauge theory. In this picture, the effect of windings of branes around is 
simulated by the rank of the gauge group. In particular, due to the Dirac quantization 
condition on fluxes, a U{n) flux on, say T^: 

F^-In, (2.19) 
n 

with In being the n-dimensional identity matrix, and {n, m) relatively prime, represents a 
single brane wound n times around with flux quantum m and resulting gauge symmetry 
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being only U (1). On the other hand, if m is an integer multiple of n such that m — pn, then 
each of the entries in the identity matrix represents a well defined U{1) flux of quantum 
p and the gauge symmetry is U{n), given by a stack of n such magnetized branes, as 
described in the last paragraph. It turns out that explicit realization of fluxes with {n, m) 
relatively prime, needs gauge configurations with non-abelian Wilson lines. 

The wavefunctions of the chiral fermion bifundamentals, with both abelian and non- 
abclian Wilson lines, involved in Yukawa computations, are given in [5] for the case of 
the factorized tori, eq. (1.1), and diagonal fluxes. For oblique fluxes, we postpone the 
discussion of non-abelian Wilson lines and rational fluxes to the last part of the paper 
(Section 7) and for the moment we consider the case of integral fluxes only. This restriction, 
nevertheless, allows for a rich structure of phenomenological value, since semi-realistic 
models with three generations of chiral fermions and stabilized moduli can be built even in 
the context of such integral fluxes, by turning on NS-NS antisymmetric tensor background. 
For example, a three generation SU{5) GUT with stabilized moduli given in [27] was 
constructed with all winding numbers, n = 1, for different stacks of branes. Also, the 
presence of a half-integral NS-NS antisymmetric tensor docs not modify any of our results, 
since all the relevant chiral fermion wavefunctions depend on the difference of fluxes along 
pairs of brane stacks which is always integral. 

3 Yukawa computation on factorized tori 
3.1 Wavef unction 

A detail discussion of the chiral fermion wavefunctions in the presence of constant gauge 
fluxes is given in Appendix A for general tori and fluxes. In the case of factorized tori, 

eq. (1.1), the six dimensional chiral/anti-chiral wavefunctions are written as a product of 
wavefunctions on T^. To show this explicitly, we present (as in Appendix A) the case of 
as an example, with case working out in a similar fashion. More precisely, considering 
that on two T^'s, fermion wavefunctions 




(3.1) 
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with their internal U (ni) x U (712) structure being represented in a manner as in eq. (A. 9), 
satisfy the equations: 







= 0, 




- A').,i.^y 


= 0, 


d2^pf + {A^ 


- 


= 0, 


d2^'L^ + {A^ 


- A%^^^^ 


= 0. 



(3.2) 

T'^ fermion wavefunctions are then constructed through a direct product of ■^^ and -0^ (in 
the notations of Appendix A): 



(1) 



(3.3) 



In particular, 



.,(2) 



(3.4) 



satisfies precisely the equations (A. 11) for chiral fermions on T^. We can further extend 
these results to show that chiral wavefunctions can also be written as a product of the 
chiral wavefunctions on three T^'s in the decomposition (1.1). 

Yukawa interaction on is then also given by an expression which is a direct product 
of the interaction terms for the three T^'s. Wavefunctions for the chiral fermions on a 
(with coordinates x, y) are expressed in terms of the basis wavefunctions tf)^'^ [5]: 



3_ 
N 





{Nz,Nt), 



j=0,...,A^-l , 



(3.5) 



with N denoting the difference of the U {ua) and U (rib) magnetic gauge fluxes as given in 
eq. (2.15), turned on along the Cartan generators, representing stacks of Ua and branes 
respectively and gives the degeneracy of the chiral fermions: 



N ^ma-mh = lab, 



(3.6) 



with nia and mj, being the 1st Chern number of fluxes along stacks a and b, with unit 
windings, as deflned through eq. (B.4). 
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Using such a basis, the chiral and anti-chiral (left and right handed fermions) basis 
wavef unctions: 



(3.7) 




are given by: 



= V'^-'^(r, ^ + c), (KT = r'^-'^if, z+o, 



= V'^-.^(f, z + O, i^ir = r^'-'^ir, z + C), (3.8) 
and satisfy the equations: 

with ( representing the Wilson hues. In the following we set the Wilson lines C = 0. 

Furthermore, expressions of the chiral and anti-chiral solutions, as given in eqs. (3-8) 
and (3.5), are well defined provided > for the wavefunctions and < for the 
wavefunctions ip-L- In these cases, for and to be properly normalized: 

/ dzdzil^iiiljiy = Sjk, (3.10) 
7r2 



= ] (3-11) 



an additional factor 

f2ImT\N\\^ 
V A" 

needs to be introduced, with A being the area of the T^. 

In fact, the basis functions (3.5) are also eigenfunctions of the Laplacian. We elaborate 
on this point more in Section 5.4 and now proceed to make use of these fermion and 
boson basis functions to determine the Yukawa interaction in the case of factorized tori 
and 'diagonal' fluxes. 
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3.2 Interaction for factorized tori 

We now summarize the basic results of [5] regarding the computations of Yukawa inter- 
actions. Such four dimensional interaction terms were obtained through a dimensional 
reduction of the D = 10, N = 1 super- Yang-Mills theory to four dimensions in the pres- 
ence of constant magnetic fluxes. The Yukawa coupling is given by 

Y,j,= [ ^fT^^]cl>U,,,, (3.12) 

JM 

where M. is the internal space on which the gauge theory has been compactified and 
and are the internal zero mode fluctuations of the gaugino and Yang-Mills fields with 
fabc being the structure constants of the higher dimensional gauge group. For the torus 
compactification that we are discussing, the internal wavefunctions are factorized into 
those depending on the coordinates of three T^'s. In turn, these involve the evaluation of 
terms of the type: 

/ dzdzTr{i)+.[(i)-,i)+\} and / dzdzTr{ij-.[(t)+,ij-]}, (3.13) 

with (f)± being the wavefunctions of the bosonic fluctuations of the ten dimensional gauge 
fields with helicity ±1 along the particular direction. Similarly ip± denotes the spinor 
fluctuations with helicities Therefore, In the factorized case of eq. (1.1), the full 

interaction term is computed as a product of three such integrals. To evaluate these 
integrals, one uses the wavefunctions (3.1) and basis functions as given in eq. (3.5). 

In the language of string construction with magnetized branes, = lab corresponds 
to the intersection number for the string starting at a stack a and ending on another one 
h. The Yukawa interaction then reads: 

Yijk = g<Jabc f dzdz'ilj''''^''{r, z).'iIj^^'^'^{t, ^).(V''''^n^, (3-14) 

with Ibc < 0, corresponding to the fact that when the intersection numbers lab and Ica 
are positive, then Ibc has to be negative, since lab + he + ha = 0. A similar expression 
exists for he > as well. To evaluate this integral, one uses an identity, satisfied by the 
theta functions appearing in the definition of the basis functions (3.5). The aim of this 
relation is to establish a connection between the wavefunctions with intersection numbers 
A^i and N2 for bifundamental states in brane intersections ab and ca with the one in the 
intersection be with = Ni + N2. However, in view of the further generalization to the 
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oblique flux case, we establish this identity explicitly in the next subsection and generalize 
it further in Section 4. 



3.3 Jacobi theta function identities 

We now explicitly prove the following theta function identity [29] used in [5] for computing 
the Yukawa couplings: 



1? 



Ni 





s 

N2 



r+s+Nim 

N1+N2 



{Zi + Z2,t{Ni + N2)) 



N2r-Nis+NiN2m 
NiN2{Ni+N2) 



{Z1N2 - Z2N1, TNiN2iNi + N2)), (3.15) 



where is the Jacobi theta-function: 



a 
b 



(3.16) 



To proceed with the proof of the above identity, we write its LHS explicitly as: 



Ni 





N2 




[Z2 



^i{^+l2?TN2 2^i{^+l2)z2 



g ^i^+'a)^ - g VJV2 



Similarly the RHS of the identity (3.15 ) can be written 



as: 



r+s+Nim 
N1+N2 





iZi + Z2,T{Ni + N2)) 



N2r-N\s+N\N2m 
NiN2{Ni+N2) 





{Z1N2 - Z2N1, tNiN2{Ni + N2)) 



mGZjv-|^_l_;V2 I3&Z UeZ 



(3.17) 



xe 



7ri 



'■(^^'^f^il^^+l^)'rN^N2(N,+N2) g2.^(^^^fj^|iL^±^ 



(3.18) 



Now, to match the zi, Z2 terms in both sides of eq. (3.15), we first note the identity: 



r + s \ f N2r — Nis \ . v f r s 

WTW2) + ''Hn,N2{n, + N2)) ^^^"^^ - '^'"'^ ^ + 7V^^^ ' ' 



(3.19) 
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and find coefficients Pi,P2, ?2 such that, 



{pih + P2I2) {zi + Z2) + {qih + 52^2) {Z1N2 - Z2N1) = {hzi + ^2^2) • 



(3.20) 



Eq. (3.20) leads to the following values for pi,P2, Qi, Q2 '■ 



Pi 



No 



P2 



ATi + iVs' N1 + N2 

1 -1 

?2 



Ni + N2' A^i + A^2 

Then the two terms, containing zi, Z2, in the RHS of eq. (3.17) can be rewritten as: 



(3.21) 



^2ni{^^+h)z^ ^2nii^+h)z2 ^ ^2niij^ + j^^+j^){z^+Z2) ^^i( ^J^^^^^ +±L_l^)iz^N,-z,m)_ 



Similarly, coefficients p, q satisfying identity: 



P 



r + s ^ Nik ^ N2I2 
Ni + N2 Ni + N2 Ni + N2 



+ q 



N2r - Nis h-h 

NiN2{Ni + N2) Ni + N2_ 
2 



Ni 



+ h 



iVi + 



No 



+ I2 



(3.22) 



= (3.23) 



No 



are given as: 



p^Ni + N2, q = NiN2{Ni + N2). 



(3.24) 



Using eqs. (3.19), (3.20), (3.22) and (3.24), the RHS of eq. (3.17) (appearing in the LHS 
of eq. (3.15) ) can be re-written : 



'f«(ir+'i)''-rAfi 27ri(j^+Zi)2i wii^+l2)^TN2 27Ti{±+l2)z2 



e "1 



EE 



EE 

Zisz hez 

g'^H -;Vi+Af2 ^ JVi+iV2 ^ Afi+-iV2 ^W^l+"2; ^^l^z^ -I- -I- j\Z\^Z2) _ 



Proving the identity, eq. (3.15), now amounts to showing that the RHS of eq. (3.18) 
matches precisely with that of eq. (3.25) with m in eq. (3.18) taking value as m = 
0,1, ,(iVi + iV2- 1). We note: 

1. When li — I2 in eq. (3.25), the terms in the RHS are identical to those in the RHS 
of eq. (3.18), with m = 0, Z4 = 0, if we identify I2 with Z3. 

When li — I2 + 1, the terms in eq. (3.25) exactly match with those in eq. (3.18) 



^^"■'\Ni+N2 ' N1+N2 ' N1+N2' 
N2r-Nis , I1-I2 
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obtained for the values m — l,l4 — with the identification of I2 with l^. 
This goes on up to h — I2 + {Ni + N2 — 1) which corresponds to the case for 
/s), m = (A^i + A^2 - 1) and k = 0. 

2. The terms obtained in cq. (3.25) for h = I2 + {N1 + N2) corresponds to m = 0, ^4 = 1 
and I2 + Ni identified with /a in eq. (3.18). 

When li — I2 + {Ni + N2) + 1 the terms correspond to the case m — l^U — 1 and 
I2 + Ni identified with ^3 in eq. (3.18). 

This goes on till h — I2 + {Ni + N2) + {Ni + N2 — 1) when they correspond to 
m = (A^i + A^2 - 1), ^4 = 1 and h + Ni identified with I3 in eq. (3.18). 

3. Similarly the terms for li — l2 + 2{Ni + N2) correspond to the terms for m — O^U — 2 
and I3 — {I2 + 2Ni) . And so on.... 

We have therefore shown a one-to-one correspondence between the terms in the RHS 
of eqs. (3.18) and (3.25). The identity eq. (3.15 ) has thus been proved explicitly. 



3.4 Application to Yukawa computation for factorized tori 



We now make use of the above Jacobi thcta identity as well as of the explicit forms of the 
fermion and scalar wavefunctions, defined in terms of the basis functions in eq. (3.5) to 
write the expression for the Yukawa interaction term. More precisely, in order to evaluate 
the Yukawa coupling given in eq. (3.14), one uses the theta identity of eq. (3.15) and the 
basis function in eq. (3.5) and proceeds by writing down: 



X 







j_ 
N2 



{N2Z,N2t), Z^0,...,N,-1, J =0,...,iV2-l. 



(3.26) 



where we have also made use of the normalization factor, jV given in eq. (3.11), and 
identified for a compactification: 



Nl = Iab, N2 = Ica, 



(3.27) 



with 



lab ^nia-mb, etc. 



(3.28) 
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giving 



N3^{N, + N2)^Icb. 

Now, using the theta identity (3.15), eq. (3.26) can be rewritten in the form: 



(3.29) 



2/mr\ 4 / lablca^ " 



Lch 



Icq i — IghJ -j-Igb^ca m 
^ab^ca^cb 





(0, rlabhal 



cb 



(3.30) 



The Yukawa interaction (3.14), is then evaluated using the orthogonahty property of the 
wavefunctions given in eq. (3.10) and reads: 



^ijk ^abcQ 



/ 2/mr y flabhaV 

V A' ) 



cb 



^ab^ca^cb 



After imposing the Kronecker delta constraint, we obtain: 



(llmry flablcaY . 



The final answer can be expressed as : 







{0,TlabIcaIcb)- 
(3.31) 



{0,TlabIcaIcb). (3.32) 



^ijk ^abc9 



l-cb 







{0,r\IabIbJca\), 



with 



_ i j k 
Oijk ^ -j— + ^ + -j-- 

^ab ^ca ^bc 



(3.33) 



(3.34) 



The result can be easily extended to the case of factorized (1.1) and the interaction is 
then written in terms of the products of theta functions of the type appearing in eq. (3.33). 
We refer the reader to [5] for the details and now go on to the generalization when fluxes 
of both oblique and diagonal forms are present. Such magnetic fluxes do not respect the 
factorization and hence involve the wavefunctions written in terms of the general Riemann 
theta functions. 



4 General tori and 'oblique' fluxes 

Let us now consider the more general case where the 2n-dimensional torus is not neces- 
sarily factorizable. A generic flat 2n-dimensional torus, T^" ~ C"/A, inherits a complex 
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structure from the covering space C". Its geometry can hence be described in terms of a 
Kahler metric and complex structure as 



ds^ 

dzi" = dxi' + r'^dy'' 



(4.1) 



where x'^,y^ E (0,1), = l,...,n, parametrize the 2n vectors of the lattice A. The 
natural generalization of the Jacobi theta function (3.16) to this higher-dimensional tori 
is known as Riemann -f^-f unctions, as defined in eq. (A. 15): 



a 
b 



i-n-{l+a) .n.{l+a) 2m{l+a) .(v+b) 



(4.2) 



As already elaborated upon earlier, in our case, although the geometry itself may be 
such that is factorizable as in eq. (1.1), the fluxes turned on, may violate in general the 
factorizable structure of the tori. Indeed, the general wavefunctions for bifundamentals 
given in terms of basis functions (A. 14)^: 



l/j^'^{z,Q) = TV . e{»^[N.i].(N./mQ)-i/m[N.^]} . ^ 



J 





(N.^IN.O), 



_ J\J- . gi7r[N-5]-(Imf2)-i-Imz _ ^ 



J 





(N-;?|N-0), 



(4.3) 



with N's being the intersection matrices, depend on such fluxes explicitly in terms of its 
modular parameter argument: Nf2; this breaks in general the factorized structure, even 
if the complex structure fl is diagonal. The explicit form of the normalization factor J\f 
appearing in eq. (4.3) is given eq. (A. 17). One needs to obtain an overlap integral of three 
basis functions of the type (4.3), in order to generalize the results of Yukawa computations 
given in eqs. (3.14), (3.30) - (3.33). 



^See Appendix A for more details on the properties of the wavefunctions and Section 2, as well as 
Appendix C for discussion on fluxes. 
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4.1 Riemann theta function identity 



We now generalize eq. (3.15) to the case of general Riemann theta functions given in eq. 
(4.2). Explicitly, we consider the LHS of our identity to be given by an expression: 



Ji 




(il|Ni -O) •'(? 



J2 





(i^2|N2-0) 



(4.4) 



where is an nxn complex matrix and Ni, N2 are nxn integer- valued symmetric matrices 
satisfying the constraints (2.17). These constraints, in turn, follow from the convergence 
of theta series expansion, as well as from the holomorphicity of fluxes: for instance, eq. 
(2.12) when p^x and Pyy components of fluxes are zero, with (i = 1,2,3) denoting 

the coordinates of three T^'s in the decomposition (1.1) and (4.1). Generalization to the 
case when p^^ and Pyy flux components are also present is discussed in Section 2, as well as 
later on in subsection 4.7, and is relevant for evaluating the Yukawa couplings in models 
with moduli stabilization, such as the one of [27]. 

Initially, we also restrict ourselves to the case when Q = rln with 7„ being a n x 
n identity matrix, implying that the geometric structure is factorized as in eq. (1.1). 
However, in Section 4.6, we generalize the results further to the case when D, is an arbitrary 
matrix satisfying the F(2,o) = supersymmetry condition, as given in eqs. (2.7) and (2.8). 
Then, using the definition of Riemann i^-functions (4.2), the expression in eq. (4.4) can 
be expanded as: 



Ji 




(illNir)-^ 



J2 





ll,l2&Z^ 



^■!Ti{j2+l2)-'N2T-{j2+h)„2m{j2+h)-Z2 



(4.5) 



Now, by defining 2n-dimensional vectors: 

d+ = 

and the 2n x 2n dimensional matrix: 




Z2 



(4.6) 



Nir 
NaT 



(4.7) 
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eq. (4.5) can be re-written as: 



3i 




(i;|NiT)-t? 



32 





[Z2 



,7ri(j+r)^ -Q- (j+f) „27ri(j+r)^ -z 



(4.8) 



Our aim in combining the terms into 2n dimensional vectors and matrices is to gen- 
eralize the procedure outlined in [29] to our situation, namely when two theta functions 
appearing in the LHS of the identity (that we propose below) carry independent modular 
parameter matrices Nir and NaT, which generally may not commute. Note that the re- 
sults of [29] are insufficient to give such an identity as they involve theta functions whose 
modular parameter matrices are proportional to each other and therefore commute. In 
order to proceed, we note that using a transformation matrix: 



T 




^ . 1 Nr^a^ 
1 -N2-'a^ 



and 



T-i = (Nr^ + N2"')-' 



(4.9) 



(4.10) 



(4.11) 



with a being an arbitrary matrix (to be determined below) and Ni,N2 being real sym- 
metric matrices, due to the condition (2.17) (for Vt = rln), one obtains: 



Q' = T Q = 



(Ni + N2)t 







a(Nr^ + N2'^)™^ 

In the following we also make use of the identities: 

(Nr^ + N2-') = Nr^(Ni + N2)N2-' = N2-'(Ni + N2)Nr^ 



(4.12) 



and 



(Nr^ + N2- 



Ni(Ni + N2)-^N2 = N2(Ni + N2)-^Ni 



(4.13) 



(4.14) 



in simplifying certain expressions. 

The transformation matrix T defined above is used to transform the product of theta 
functions in the LHS of eq. (4.8), in terms of a finite sum over another product of theta's. 
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now with modular parameter matrices: (Ni +N2)r and q;(Ni ^ + N2 ^jra^. Explicitly, 
we can write the terms appearing in the exponents in the RHS of eq. (4.8) as: 

(J + If-Q- (J+ I) = (J+ If ■ (T-'T) ■ Q • {T^{T-'f) ■ (J+ I) (4.15) 
(J+ f)^ . z = (J+ l)^(T-iT) . z. (4.16) 



Then using: 



T-z 



Zl + Z2 



(4.17) 



and 



(ji + /i)(Ni-^ + N2-^)-^N2-^ + (j2 + /2)(Nr^ + N2-^)-iNl 



1 _L ivr-^-iivr-1 



((ji + ^i) - O2 + /2))(Nr^ + Nj^)-ia-i 



,(4.18) 



we can re-write eq. (4.5) as. 



(q;-1)^(Ni + N2 ')-M(ji + ^i) - (J2 + I2)] 



Jl 




32 





(i^2|A^2T) 



g7ri[{((/i+ii)Ni+02+i2)N2)(Ni+N2)-i}-(Ni+N2)r-{(Ni+N2)-HNi(/i+ii)+N202+i2))}] 

/I,Z2eZ" 

^g27ri{[((/i+/i)Ni+02+i2)N2)(Ni+N2)-^]-[z-i+z-2]} ^ 
^7ri{[(0i-j2)+(i7-/2))Ni(Ni+N2)-iN2a-i]-[a(Ni-i(Ni+N2)N2-i)Tar].[(a-i)TN2(Ni+N2)-iNi(0i-^^ 



Xe2'''{[((-'^i--?2)+{i7-/2))Ni(Ni+N2)-iN2a-i]-[aArj-iil-aiV2-iz2]} 



(4.20) 



Now, to reexpress the above series expansion in terms of a sum over theta functions 
with modular parameter matrices: Ni + N2 and q;(Ni~'^ + N2~^)q;^, we rearrange the 

— * — * 

series in eq. (4.20) in terms of new summation variables ^3, ^4, m, whose values and ranges 
will be assigned later. In the course of going from eq. (4.20) to (4.22) below, however, one 
needs to make sure that such redefined variables are integers. This requirement constrains 
the matrix a whose 'minimal' solution will be taken to be 



a = (detNi det N2)/. 



(4.21) 
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We will later on discuss also the possibility of choosing other forms of a and show that 
such choices lead to the cyclicity of the superpotential coefficients, as in eqs. (3.33), (3.34). 
Using eq. (4.21), the RHS of eq. (4.20) takes the form: 

^g7ri[(/iNi+j2N2+mNi)(Ni+N2)-i+Z3]-(Ni+N2)T-[(Ni+N2)-HNi/i+N2j2+Nim^ 
„27ri[0iNi+j2N2+mNi)(Ni+N2)-l+Z3]-[2l+22] y 



dct N]^ dot N2 



,27ri[0-i-j2+m) '^^det det N^^ +^4]-det Ni det NalNr^il-N-^^a] 



(4.22) 



This series can now be reexpressed in terms of a finite sum over product of generalized 
theta functions given in eq. (4.2), leading to a generalization of the identity (3.15) to: 







(^1|Nit)-^ 



J2 





(i^lNar) 



(/iNi + /2N2 + m.Ni) (Ni + Na)-^ 


[(ji - J2) + m] — ^ ^ 



(iI + £2|(Ni + N2)r)x 



det Ni det N2 







((detNidet N2)(Ni-^^l -N2-^i^2)|(detNidetN2)^(Nr^(Ni + N2)N2-^)r), 

(4.23) 

where fh = miel are all vectors generated by the basis vectors el: 



etc., 



(4.24) 



and lied within the unit-cell defined by the new basis vectors: 

e' = e(detNi det N2)(Nri(Ni + N2)N2-'). 



(4.25) 



The above identity already assumes the form Q = t/„ for the complex structure of 
T^". As mentioned already, in subsection 4.6 below, we make further generahzation to 
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include arbitrary complex structure fl as well. Also, note that, due to the identities (4.13) 
and (4.14), the theta functions appearing in the RHS of eq. (4.23) satisfy the constraint 
(A. 16) with respect to their own arguments. 

4.2 Proof of the identity 

We now show the equality of the series expansions (4.20) and (4.22) to estabhsh the 
identity eq. (4.23). We also show that matrix a needs to be chosen as in eq. (4.21) for 
showing the equality of eqs. (4.20) and (4.22) for the case when detNi and detN2 are 
relatively prime. In other cases a can be chosen as the least common multiple of det Ni 
and detN2. Here we assume them to be relatively prime, while the remaining cases can 
be worked out in a similar fashion. 

We now follow an exercise similar to the one in Section 3.3, to show that series in eqs. 
(4.20) and (4.22) precisely match with m restricted to be an integer, provided a is given 
byeq. (4.21). 

1. When /i = I2 in eq. (4.20), we have: 

(ZiNi + Z2N2)(Ni + N2)-' = /2 (4.26) 

and 

{h - /2)Ni(Ni + N2)-'N2Q;-^ = (4.27) 

These terms are exactly same if we consider the series given in eq. (4.22) for the 
values /3(= ^2), ^4 = and m = 0, irrespective of the choice for the matrix a. 

2. In order to see the restriction on the matrix a, one needs to understand how the 
nonzero integers /4 7^ in eq. (4.22) are generated from the terms in eq. (4.20). In 
other words, one needs to make sure that 

(h - i2)Ni(Ni + N2)-^N2Q;-^ = l1 (4.28) 

is an integer. This in turn is possible only if Zi is of the form: 

/i = /2 + /4«N2i(Ni + N2)N^^ (4.29) 

However, since I4, Ni, N2, take integral values, the RHS in eq. (4.29) is an integer 
only if ci;(N]"^ + ^) is an integer. In other words, for det Ni and detN2 relatively 
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prime, a needs to be of the form: 



a = (det Ni det NajP. (4.30) 

with P being an arbitrary invertible integer matrix. 'Minimal' choice also demands 
detP = 1, otherwise Z4 will not span over all integers. Then, since P is invertible, it 
is fixed to be the identity matrix. We have therefore established the restriction on a 
as in eq. (4.21). At the same time, we have also proved that the series in eqs. (4.20) 
and (4.22) precisely match for m = provided I2 + detNi detN2i4N^^ is identified 
with I3 in eq.(4.22). Note that (det Ni det N2)N2 ^ is also integer valued and ensures 

— * 

that such an identification with I3 holds. 
3. On the other hand. When li — I2 + m in eq. (4.20), we end up with terms like: 

{hN^ + /2N2)(Ni + = l2 + m.Ni(Ni + Na)"^ (4.31) 

and 

~ detNidetN2 " detN^detN^ ^^'^^^ 
These terms can also be obtained in the series (4.22), for the following values of the 

— * — * — * 

variables: l3{= I2), Z4 = 0, m arbitrary. However, as seen above in eqs. (4.28), (4.29), 
the sum over m is finite due to the fact that 

LdetNidetN2N2^(Ni + N2)Nr\ (4.33) 

for L arbitrary integers, contributes to the values of U in the RHS of eq. (4.22) by 

— * — * — * ~* 1 

an amount L, while setting m to zero. Is is identified with I2 + detNi detN2LN^ . 
In other words, we have shown that the sum over m in (4.22) is over all integrally 
defined vectors in the unit cell generated by the basis elements: 

e' = e det Ni det N2N^i(Ni + N2)N^-^ (4.34) 

with e being the elements of the canonical basis (4.24). 

We have therefore proved that identity eq. (4.23) holds by explicitly showing a one to 
one correspondence between the series in eqs. (4.20) and (4.22). 
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4.3 Yukawa expressions for oblique fluxes 



We now use the wavefunctions given in eqs. (4.3) and (4.2), to obtain the expression of 
Yukawa interactions when obhque fluxes, specified by intersection matrices 



Ni = F„ - Fb, N2 = Fe - Fa, N3 = Fe - F,. 



(4.35) 



are turned on along branes a, b and c. As aheady mentioned, in eq. (2.15), Ni, N2 
and N3 are all real symmetric matrices (in the absence of components p^x, Pyy) and in 
addition the complex structure matrix is chosen to be proportional to the identity: t/„, 
with T complex. We then have: 



/'^^(z, n = tQ ■ i/j^'^^iz, n = Tin) = (2?) (yoZ(T2")) ' (I det Ni|. I detNa I (/mr)^) 



^ giTrNa.^Im z/Im 







(Ni • f iNi • r) ■ ?9 



J 




(N2-f|N2-T). (4.36) 



Using the Riemann theta identity derived earher in eq. (4.23), eq. (4.36) can be rewritten 
as: 



m 

^(iNi+jN2+?nNi).N3-i,N3^^^ . ^ 



(IdetNil 


|dct N2I (/mr)^) 




det N3 





X 



[{i - 3) + m] 




N1N3-1N2 

det Ni det N2 



(0|(detNidetN2)'(Nr^N3N2-')T). (4.37) 

Note that the integrality condition (A.16) is maintained by -0(?Ni+jN2+mNi)N3-i,N3(£j 
appearing in the RHS of the above equation, since the expression 



(iNi+jN2 + mNi)N3-' 



(4.38) 



is always an integer. On the other hand, the sum m in eq. (4.37) is over the integers 
inside the cell generated by the lattice vectors in eq. (4.34) and total number of them is 
given by the volume of this compact space. The size of the cell, i.e., its volume matches 
with those in eq. (3.29) and (3.30) for the case which is just the number, A^3 = I^b in 
eq. (3.29), of chiral states for brane intersection be. However, the situation is different for 
T^", n > 1. This becomes clear by observing that the size of the cell given in eq. (4.34) 

— * 

is bigger than the number of states {k) in the intersection N3 between the branes b and 
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c by a factor ciei(detNi detN2N2 ^N^ This factor, on the other hand, for is unity. 
We therefore notice that the sum m is over many more terms^ than the actual number of 

— * 

states (k) in the intersection N3 between the branes b and c. 

The extra factor of terms appearing in eq. (4.37) can be explained by noticing that the 
sum over terms in eqs. (4.37) and (4.39) is over the states ^(^Ni+/N2+mNi).N3-i,N3^^) ^^^^^ 
are inside the cell in eq. (4.34) and contribute to the Yukawa coupling by the orthogonality 
relation eq. (A. 18). As any state (with more details given in the subsection-4.4) k, 
satisfying integrality conditions such as (A. 16) is defined only upto the integer lattice 
shifts, one therefore has appearance of the same states inside the volume of lattice (4.34), 
multiple times. In other words, for any given state, in the RHS of eqs. (4.37), all those 
integer vector (m) shifts also contribute to the sum which satisfy the integrality condition 
for mNiNa"^ inside the cell (4.34). Exphcit solution of this condition is presented later 
on in section 4.4 in eq. (4.45). 

Then, as in the case, orthonormality of wavefunctions (A. 18), implies that the 
Yukawa coupling, whose explicit form is given in section 4.4, can be 'formally' written in 
a form : 



(IdetNil 


|detN2|(/mT)3) 




detNg 





^ '^fc,N3-i(Nii+N2j+Nim) 



det Ni det N2 





me{e'} 

(0|(detNidetN2)'(Nr^N3N2-')T), (4.39) 



where by the summation index rh e {e'}, one means to sum over all integer points inside 
the lattice generated by e'l, e^ - ■ ■ in eq. (4.34) and the Kronecker delta is to identify all 
the states k upto integer shifts. 

The above expression reduces in the case of flux compactification to eq. (3.32), since 
the Kronecker delta constraint has a unique solution in such a situation. To compare the 
two expressions, note that the indices i, j. A; in the factorized case are scaled with respect 
to the one of general tori, by the factors and respectively. Then, the Kronecker 

delta constraint in eq. (4.39) precisely matches with the one in eq. (3.31). In the case of 
general tori, however, the constraint implies that the interaction terms involve the states 



^We thank R. Russo and S. Sciuto for invaluable suggestions on this point. 
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which satisfy the equation 



Ng/c = (Nii + Nai + Nim) 



(4.40) 



among the vectors N^i, Naj, N3A; for m integers inside the unit cell given in eq. (4.25) 

— * 

and corresponding states k are only defined upto integer lattice shifts. We now find all 



such solutions of the lattice shifts in the next subsection and present the explicit answer 
for the Yukawa coupling for general tori. 

4.4 Explict Yukawa coupling expressions 

In this subsection we now present the set of terms that contribute to eqs. (4.37) and 
(4.39). In order to clarify the situation we analyze the correspondence between the chiral 
multiplet families of states such as the ones appearing in eq. (4.38) and the fluxes along 
the branes. Our discussion is restricted to N being real symmetric matrices, due to the 
imposition of the Riemann conditions (2.17) for the special complex structure Q = r/„ 
under discussion. 

For a given pair of brane-stacks with intersection matrix N, the condition eq. (A. 16) 
that a state i needs to satisfy is N.i — integer. The solution of this condition is: i — N~^e, 
with e being the integer basis vectors in an n- dimensional space as given in eq. (4.24). 
The states are therefore generated by the set of n vectors: ii = e^N"^, with subscript 
i = 1,2 ■ ■ ■ n and are det(N) in number, namely those which are inside the cell generated 
by e^'s. Here and in following we also keep in mind that all the chiral multiplet states that 
we are discussing, are defined only upto the shift by integer lattice vectors e^'s. 

To give an example: for n = 2 (corresponding to T^), with 




(4.41) 



we have the basis vectors for generating the states: 



ii = 



(a/3 -7^) 



1 




) 



?2 = 



{aP - 7^) 



1 




(4.42) 



To obtain the degeneracy count, we note that for the above example we have: 



ei = ail + 7^2, 
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62 = 7^1 +/5«2- 



(4.43) 



The number of independent states inside the cell with lattice vectors el and 62 is then the 
determinant of the above transformation which is detN. A generic state appearing in eq. 
(4.40) then has a form: 

i = mill + 7712^2, 3 = niji + 712^2, k = pih + p2/c2- (4.44) 

with j'j, ki defined in a similar way as in eq. (4.42) with respect to the corresponding 
intersection matrices. Also, integers mi,ni,Pi label the states of a chiral family in a given 
brane stack. 

We now go on to give explicit solution for the vector fh that contribute to the sum 
of terms in Yukawa coupling expressions (4.37) and (4.39), namely those inside the cell 
defined in eq. (4.34). The size of the cell, namely the number of states that it contains 
is equal to det(detNi detN2N2 ^(Ni +N2)N]^^), as stated earlier. In a situation with 
2x2 matrices, for example, it is det'N-idetN2detN3. For illustration purposes we restrict 
ourselves to the discussion with 2x2 matrices. However, all the results we write below 
are valid for other situations as well. 

Now, restricting to this 2x2 case for the simplicity of discussion, we write all possible 
solutions for m that provide integer solutions for mNiNs^^, as appearing in the definition 
of states in eqs. (4.37), (4.38), and show that they are detN-idetN2 in number. So that 
the degeneracy of the state matches with detNidetN2detNs given in the last paragraph. 
To compare, note that for a diagonal flux situation, as in section-3, we have m — as 
a single solution of an analogous condition mriin^^ — integer, corresponding to the state 
degeneracy which is n^. 

The integer solutions for rnNiNa"^ are: 

m = pHetNiNgNr^ + pdetN2NsN2''\ (4.45) 

where p is all integer vectors within a cell generated by edetJ>i2^2~^ and p is all inte- 
ger vectors within a cell generated by edetN It is easy to see that m satisfies 
777.NiN3~^ = integer (by making use of Ni = N3 — N2). Together, for every solution of 
the first term in m we have detNi solution for the second term and this goes on for detN2 
number of terms from the first term. So that total degeneracy of such m is detNidetN2, 
as stated earher. 
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About the states: m given in eq. (4.45) defines a periodic set, in the same way as for 
the case m — ris defines the periodic set of states in the RHS of eqs. (3.30) and (3.31). 
There the states are exphcitly given as A; = (0), {ni/n^), (2711/713), • • • [(713 — 1)711/713] with 
a periodicity 713 for this series. Various states inside the cell (4.34) can also be found using 
eq. (4.40) and making use of the condition: Ni = N3 — N2 as: (also the fact that any 
state is defined upto integer vectors). The states are: 



k = pdetNiNs'^ + pdetN2Ns'''^ etc. 



(4.46) 



and the state degeneracy is detNidetN2detNs. 

The Yukawa coupling can now be written in an explicit form given by a sum of 
detN-idetN2 number of terms, which can be read off from eq. (4.37) directly, with m 
replaced by 

m + pHeiNiNgNr^ + pdeiNaNgNa"^ (4.47) 

and now such m are the unique solutions of eq. (4.40) where all other solutions defined 
upto the shifts in m by pdetNiN^N^^ + pdetN2N^N2^ have been identified. 
Eq. (4.39) now reads as: 



Yijk = 9<Ja,c (2t)^(yo/(T2"))- 



(IdetNil 


|detN2|(/7nT)^) 




detNg 





[{(?- J) + (^Ns - ?Ni - jN2)Ni-};™g- + (P^ +#^)] 



pS 

N2 



(0|(detNidetN2)'(Nr^N3N2-V), 



(4.48) 



or equivalently: 



•y, 
ijk 



gCTabc 



(2t)^ (\/o/(T2«))" 



(IdetNil 


|detN2|(/7nT)3) 




det N3 





p,p 



i-J + k) 



N2 



detNidetN2 







(0|(detNidetN2)^(NriN3N2"V). (4.49) 



Note that the sum over m is now broken into sum over p and p. We end this discussion 
by reminding ourselves once again that p runs over all the states inside the cell generated 
by eidetN2^2~^ and e^deiN2N2~'^. Similarly p runs over all the states inside the cell 
generated by el(ietNiNi~^ and e^detNiNi"^. 
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We now present two explicit examples, one for the oblique situation and the other for 
the commuting diagonal fluxes. We show that our answer for the diagonal flux is identical 
to the one for the diagonal yukawa couphng expression given in [5] for T^". In fact this 
holds for any set of fluxes with Ni, N2, N3 diagonal. On the other hand, we also show 
that the set of terms given above in eqs. (4.48) and (4.49) can also be summed up in a 
number of cases, for the oblique cases as well. 



Example : Oblique flux 



For the oblique case, by taking two noncommuting matrices: 

Ni = 




(4.50) 



we have: 



(deiNi)Ni-^ = 




, (defN2)N2-' = 



(4.51) 



The set of integer points inside the cell generated by el(ietN2N2"^ = (2, 0) and ei(ietN2N2~^ 
(0, 1), are: (0,0) and (1,0), as det{det]Si2^2^^) = 2. The set of integer points inside the 
cell generated by eldetNiNr^ = (2,-1) and e^detNiNi"^ = (-1,2), are : (0,0), (1,0) 
and (0, 1), as (iet((ietNiNr^) = 3.^ 

Now, to illustrate our method, we concentrate on finding a particular Yukawa interac- 
tion among states: i — j — k — (0,0). This particular Yukawa now has the form, making 
use of Eq. (4.48) as: 



1^000 = 5<7ooo {2^y iVol(T'^)y 



(IdetNil 


|detN2|(/mr)2) 




det N3 





^detN2 "'"^detNi)] 





(0| (det Ni det N2)'(Nr^N3N2- V)), 



^Another example with mixed eigenvalues for the matrix Ni can be constructed by exchanging the 
ofF-diagonal and diagonal entries in eq. (4.50) for Ni. Such an example will be relevant for the situtation 
discussed in later sections where intersection matrices with both positive and negative eigenvalues are 
discussed. 
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To see what terms in p and p dependent arguements appear in the sum, we write down 
all the possibilities that arise from the combinations: 



[p + p- 



detN2 detNi ' 



p 



-,12 1 
1 3 I 1 2 



(4.52) 



with p = (0, 0), (1, 0) and p = (0, 0), (0, 1), (1, 0). All the six possibilities then imply that 
in Theta function we get the following explicit sum: 



+ 





' [(0,0)] " 








'[(hi)]' 




'[(hi)]' 
































1? 



nil)] 







(0| (det Ni det N2)^(Ni-^N3N2- V)) 



(4.53) 



where a common modular parameter arguement of the all the six Theta terms have been 
written outside of the bracket for saving space. The integer sums of the six terms over 
integer / are of the forms: 



y^g[r+(gi,g2)](detNidetN2)^(Ni-iN3N2-V)[r+(qi,g2)] 

r 

with I + {qi,q2) given explicitly as: 

r+(o,o), r+(^,o), r+(^,^), r+(^,^), r+(^,^), r+(^,^), 

for the six terms in eq. (4.53). It can also be seen that we can write them as: 



(4.54) 



(4.55) 



h 



+ 



m I 2n 
2 3 

n 
3 



li\ 1 / 3 4 \ / m 
lol 6 1 2 / \ n 



(4.56) 



with m — 0,1 and n — 0,1, 2. Now, using the inverse of the matrix 



1 /3 4 



6 \0 2 



(4.57) 



appearing in eq. (4.56): 



P-1 = 



2 -4 
3 



(4.58) 
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we can write eq. (4.56) as: 




(4.59) 

I 2y y '6I2 J \^ J 

with m — 0,1 and n — 0,1, 2. 

It can now be seen that as /i, I2 vary over all integers, and m = 0,1 and n = 0, 1, 2, 
then the combination of terms in the big square bracket in eq. (4.59) also span over ALL 
integers. As a result we are able to take the factor of matrix P out by summing over all 
the six terms, while reducing the six terms in eq.(4.53) to one. The net result is then the 
arguement of theta function modifies by the factor: 

(det Ni det N2)^(Nr^N3N2" V) ^ P^(det Ni det N2)2(Nr^N3N2" V)P (4.60) 

and final answer for Yukawa coupling is: 

>ooo = » (2t)^(yo/(T2"))- 



(1 det Nil 


|detN2|(/mT)3) 




detNg 





X 








|P^(det Ni det N2)2(Nr^N3N2"V)P). 



We can similarly take care of other nonzero values i,j,k etc. as well, but details are being 
left. 



Example : Diagonal Flux 

We take another example, now with diagonal fluxes : 



Then: 




(4.61) 



(4.62) 



Set of integer points inside the cell generated by eldeiN2N2 ^ = (2, 0) and e^deiN2N2 — 
(0,5), are: (0,0), (0,1), (0,2), (0,3), (0,4), (1,0), (1,1), (1,2), (1,3), (1,4), as detidetN^N^'') = 
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10. On the other hand, set of integer points inside the cell generated by eldetNiNi ^ — 
(3, 0) and e^deiNiNr^ = (0, 2), are: (0, 0), (1, 0) (0, 1), (1, 1), (2, 0), (2, 1), as det{detNiNi-^) 
6. 



We now have: 



l + (p 



detN2 

which can also be written as 

N2 



+ p- 



/i 



detNi ' 



l + p 




(4.63) 



l + {p 



detN-2 



+ p- 



detNi ' 




(4.64) 



with pi = 0, 1, vi = 0, 1, 2, 3, 4, pi = 0, 1, 2, p2 = 0, 1. 

By taking a factor of (jetN^detN2 above equation can also be rewritten as: 



N1N2 



detN-idetN2 



No ^ Ni 



detN-2 



detNi ■ 





(4.65) 



with pi — 0, 1, P2 — 0, 1, 2, 3, 4, pi — 0, 1, 2, p2 = 0, 1. It can again be checked explicitly 
that it leads to ALL integer variables inside the square bracket. The net result of summing 
over different terms in the diagonal case therefore is the appearance of the matrix outside 
the square bracket: det^T^etN^' When multiplying the modular parameter arguement as 
appearing in eq. (4.48), from both left and the right, this precisely reproduces a modified 
modular parameter which matches with the known diagonal flux solultion for Yukawa 
coupling in [5]. This holds for the diagonal flux in general, not restricted to the example 
above. 



4.5 arbitrary-a 

The results, obtained so far in this section, are derived for a particular choice of a given 
in the eq. (4.21). However, all the results can be re-derived for arbitrary a, appearing 
in eq. (4.9) etc.. For the factorized case, we saw in that the Yukawa coupling expression 
(3.32) can be recast into a symmetric form in eq. (3.33) (apart from the prefactor), 
where the arguments of the Jacobi theta functions are invariant under a cyclic change: 
a — > 6 — > c. This is due to the cyclic property of the superpotential coefficients obtained by 
a third derivative of the superpotential Wijk- The prefactor does not obey in general this 
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symmetry, since it depends on the wave function normalizations (Kahler metric). Here, 
we show a similar cyclic property in the non-factorized case, given above in the Yukawa 
coupling expression (4.49), by making different choices of the matrix a. in eq. (4.21). Note 
that different choices of this matrix provide equivalent expressions for the wavef unctions, 
and in turn Yukawa couplings, since they are related though a change of variables inside 
the theta sum. The a matrix can be chosen appropriately so that the redefined variables 
in eqs. (4.29) and (4.33) are well defined integers. Below we present a few examples with 
different choices of «, to demonstrate the cyclicity mentioned above. 
Eq. (4.22), for arbitrary a, can be written as: 



i[(ilNi+i;N2+mNi)(Ni+N2)-i+i3]-(Ni+N2)r-[(Ni+N2)-HNiil+N2i;+Nim)+i3] 



xe 



27r Ni +J2N2 +mNi ) (Ni +N2 ) " ^ +^3] • [il +22] 



) 



X ( g7ri[(j^-i;+m)Ni(Ni+N2)-iN2a-i+i4]-[aN^i(Ni+N2)N2 V]a'r.[(a-i)'rN2(Ni+N2)-iNi(iI-i^ 
y 27ri[(il-i^+m)Ni(Ni+N2)-iN2a-i+i4]-[o!Nj;i2l-N2 1^2] ^ 

^ ' (4.66) 

provided ^4, defined in eq. (4.28), is an integer vector, and so is m given in eq. (4.33). In 
addition the unit-cell, within which m lie, is now defined by the basis vectors : 



d = ea(Ni-i(Ni + N2)N2-^). 
Moreover, eq. (4.23) takes the form: 



(4.67) 



1? 



Jl 




J2 





(ii|N2r) 



(il + i^2|(Ni + N2)T) 



(jiNi + J2N2 + m.Ni)(Ni + N2)-i 


[(/i - /2) + m]Ni(Ni + N2)-iN2a-i 


{a{-Nr^z-[ - N2-'r2)|a;(Nr^(Ni + N2)N2-V)q;^) . 



(4.68) 



It is then easy to see, all equations from (4.36) to (4.39) go through for arbitrary a, giving 
the following expression for the Yukawa couplings: 



(IdetNil 


|det N2I (/mr)3) 




detNg 





X 
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(0|a(NriN3N2-V)a^). 



(4.69) 



where the sum m is now over all the integer solutions of mNiNa"^ in the cell given in eq. 
(4.67). Explicit contributions to this sum, of course, will depend on the exact form of a. 
In subsection 4.4, we have presented the case of a = detN-idetN2. 

We now study how the above expression (4.69) reduces for another choice of a, such 

as: 

« = N3-^Ni(detN2.detN3). (4.70) 

Note, for this choice of a, that the degeneracy of states in the cell given in eq. (4.67) is 
det{detNsdet'N 2N2~^)- As a result, for the case of 2 x 2 matrices for example, one now 
expects the sum over m to run over detN2detNs values. Explicit solutions are now given 
as: 

m = pHetNaNgNa"^ + pdetNs, (4.71) 

where p is all integer vectors within a cell generated by edetN^Ns'^ and p is all integer 
vectors within a cell generated by edetN 2T^2~^- It is again easy to see that rh satisfies 
mNiNa"^ = integer (by making use of Ni = N3 — N2). 

The characteristic of the function in eq. (4.69), becomes: 

Nr^Ns 



(-j + k)'N2a-^ = (-A;Ni + iNi + mNi; 

(-k + i)N3 



(det Na.det N3) 



(det Ns.det N3) ' 

where in the first equality we have made use of eq. (4.40). Also we have. 



(4.72) 



a(Nr^N3N2"V)a'^ = (N3"^Ni)(Nr^N3N2"^)(NiN3"^)r(det N2. det N3)^ 

= (N2"^NiN3-V)(detN2.detN3)2. (4.73) 

The Yukawa couplings then read (following the exercise performed in subsection 4.4): 

■(|detNi|.|detN2|(/mT)3) 



p,p 



detNadetNa ^ V^^detNg ^^detNa- 



X 



detN3| 

(0|(detN2detN3)'(N2-^NiN3-^)T), 

(4.74) 
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where the summation over indices p and p is explained earher after eq. (4.71). We can 
also explicitly obtain the sums, as done for various examples in the last subsection. 

Now, a comparison of eqs. (4.49) and (4.74) shows a symmetry between the ■j?-function 

— * 

characteristics in these cases, including the summation variables p and p. It is obvious 
that the replacement i ^ j,j ^ k,k ^ i and Ni — > N2,N2 — > N3,N3 — > Ni in eq. 
(4.49) results eq. (4.74). We have thus established that just as in the factorized case, 
for oblique fluxes too, one can show the cyclicity property of the Yukawa superpotential 
coefficients, as naively expected. 

4.6 General complex structure 

In the previous subsections 4.1 - 4.3, we have confined ourselves to the complex structure 
matrix Q = rln for a 2n dimensional torus. This implies the restriction to orthogonal tori, 
a solution which is already used in many phenomenologically interesting models. However, 
the results are easily generalized to complex structure with arbitrary fl. More precisely, 
to write down an identity generalizing eq. (4.23) one starts with the product expression 
given in eq. (4.4) and rescales Ni, N2 in eqs. (4.7) - (4.25) to NiQ/r, N2Q/T. At the 
same time, the matrix a in eq. (4.21) is also rescaled : 

aQ 

a ^ a = detNidetN20/T = — . (4.75) 

T 

Moreover, one needs to take into account that in relations such as (4.10) earlier, we have 
made use of the property N-^ = N, which is true for the complex structure of the form: 
rln- Replacements: Nr — > Nil are, however, to be done in the original expression. 

Explicitly, under the changes mentioned, the transformation matrix T in eq. (4.9) 
remains unchanged, while its transposition in eq. (4.10) is now written as: 



a 



Also, (4.11) is unchanged, whereas Q' in eq. (4.12) goes over to 

^ / (Ni + N2)l^ , , , 

Q' = T • Q • = ^ ^ , 4.77 

\^ «(Ni-i + N2-^)llV ' 

where we have made use of the fact that both (Ni + N2)J1 and (Ni~^ + N2~^)il^ are 
symmetric matrices, due to the condition (2.12), with N defined in eq. (2.15). Then 
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expressions (4.17) and (4.18) remain unchanged, while (4.19) is modified to: 



+ h) + Nr^^(Nr^'^+N2-^^)-H/2 + k)] 



The identity (4.23) then takes the form: 



(4.78) 



1? 



Ji 
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(4.79) 



(/iNi + /2N2 + m.Ni)(Ni + N2)-i 


^\ I ^1 Ni(Ni+N2)~iN2 

[{31 - 32) + m\ ^ 



(il + iiKNi+Na)^) X 



det Ni det N2 







((det Ni det N2)(Nr^il - Na'^i^s) | (det Ni det N2)2(Nr^(Ni + N2)N2-'f^^)), 
leading to the expression for the Yukawa interaction: 



Yiik^oa^cg (2t)^(Vo/(T2")) 

—3 + ^) det, Ni det No + 



(IdetNil 


1 det N2 


|detn|) 




detNal 





X 



E 



detNidetN2 ' V/'detN2 ^■'^detNi' 







(0| (det Ni det N2)^(Nr^N3N2"^0^)). 

(4.80) 

We leave the rest of the details, which readers can work out. 
4.7 Hermitian intersection matrices 

In subsections 4.1, 4.2, 4.3, we have assumed that intersection matrices Ni,N2 etc. are 
real symmetric. As explained, this restriction originates from the case when fluxes 'p^xi Vyy 
are zero and the intersection matrix N is represented by the real matrix "p^y in eq. (2.15), 
which is symmetric whenever the complex structure is of the canonical form: i7 = il^- 
Moreover, the Yukawa coupling expression was generalized nicely in the last subsection to 
the case of arbitrary complex structure, as well. 

In this subsection we discuss the case when fluxes p^jj; and 'Pyy are also present, in 
addition to those of the type 'p^y and Pyx- Furthermore, all these fluxes are constrained by 
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the conditions (2.7) and (2.8) giving a resulting (1, 1) - form flux which can be represented 
by the Hermitian matrix (2.9), (2.10). We exphcitly present the case of Q = ild solution 
{Id '■ d-dimensional Identity matrix), which is particularly simple, since in this case due 
to constraints (B.l), the Hermitian flux has the simple flnal form of eq. (B.2). The 
generalization to arbitrary complex structure can also be done, but is left as an exercise 
to the reader. 

Wavefunctions on T^, as given in eq. (4.3), satisfy the following field equations (A. 12) 
and (A. 13): ^ 

dat + (A' - A%xt = 0, {% = 1, 2, 3). (4.81) 

We now show that the solution for the above equation, together with proper periodicity 
requirements on T^, is given by the basis elements: 

^^^^{^ = Myj{z,z)-e{z,z) 

(Nr • ^ |Nr • ih) (4.82) 

where Nr, is a real, symmetric matrix. 

The wavefunction given in eq. (4.82) satisfies the Dirac equations (4.81) for the fol- 
lowing gauge potentials. 

{A^ - A^),^ = (^) ^,(Nr - iN,),-, (4.83) 

which exactly matches with eq. (A. 13) for the complex structure D, — il^. The intersection 
matrix is therefore given by : 

N = Nr - iNi, (4.84) 

where we identify, 

= Ply -Ply, ^i=p:.-pL- (4.85) 

The wavefunction described in eq. (4.82) can be re-written in terms of the real coordinates 
X and y as well as matrices Nr, Nj. By a slight abuse of notation, below, only for 
this subsection, we use Nr = p^y, Ni = p^x, by setting p'^'s to zero in eq. (4.85) and 
suppressing the superscript a in p". Such a notational change, helps to make comparison 
^See Appendix A for details. 
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of the transformation rules we derive for the wavefunction written above in eq. (4.82) with 
general transition functions, consistent with the gauge transformations along the 2n non- 
contractible cycles of T", given in [5]. These transition functions are written in equations 
(4.40), (4.41) of [5] for the fields that transform in fundamental representation rather than 
as bifundamentals. Hence, the notation changes above are meant to make the expressions 
consistent with the ones of [5] . 

The wavefunction (4.82), in the real coordinates x and y, then reads: 



QMi)[{h+ji)-p^iyjiij+jj)]^2i7T[{ii+ji)-p^iyj-^ (4.86) 

This expression in terms of real coordinates is useful in comparing the transformation 
properties of the wavefunction over with the one in [5]. The transformation properties, 
as derived from eq. (4.82), are given by, 

(4.87) 

provided that 



• (NR)jj = Pxiyj e Z, i.e Nr is integrally quantized. 



• j satisfies j ■ Nr, e Z". 

We therefore notice that the integer quantization is imposed only on the symmetric part 
Nr of the intersection matrix from the periodicity of the wavefunction as well. However, 
Dirac quantization already imposes both p^y and p^x to be integral for unit windings, as 
discussed in Section 2. 

Using eq. (4.86), the expressions (4.87) can be re- written in terms of real coordinates 

as: 

ip^'^{x + n + iy} = e^''["*(P-*«^~'^-'-^^^'l • ip^'^{x + iy), (4.88) 

^^^^{x + i[y + n]) = e-'^'^'^'^-^y'-'P-^-'^'''^ -i^^'^ix + iy). (4.89) 

In order to see that eqs. (4.88) and (4.89) are the proper transformation properties of the 
fermion wavefunction over T^, let us compare them with the the transition functions eq. 
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(4.41) of [5] given for a fundamental representation in six real coordinates Xi, I — 1, - ■ ■ ,6, 
as used in our eq. (2.1) as well. After changing variables first to the coordinates 
i — 1, 2, 3 and then making coordinate transformation to described in Section 2, 

the general transition function is given by, 

-^[Xi.yi) — e*'^[("»i+«"0--'^ij(2/-''W)+(imi+ni).F-(a:-'+i2/J)] (4.90) 

In correspondence to the transformation along the 1-cycles, the integer parameters on Xi 
and Hi are denoted as mi and rii respectively. One then has two cases: 
Case -I : When ni — 0, i.e x — > {x-\-m), eq. (4.90) reduces to 

— g'^Mm-Fij-y^) ^ (4.91) 

where we used the hermiticity property of F. Using the expression (B.2 ) in eq. (4.91), 

we recover the transformation given in eq. ( 4.88). 

Case -II : When rrii — i.e y — > (y + n), eq. (4.90) takes the form. 

Again, using eq. (B.2 ) in eq. (4.92), we reproduce the transformation (4.89). 

It can also be easily seen that the basis wavefunctions given in eqs. (4.82) and (4.86) 
satisfy the orthonormality condition 

/ (^^>)tV'^^'^ = 5.,-, (4.93) 

by fixing the normalization constant to 

jVj = (2"|detNR|)'/' • Vo1(T2-)-V2^ Vj . (4.94) 

We have therefore confirmed that the wavefunction written in (4.82) is not only a solution 

of the field equation, but also has the correct periodicity properties on the torus under 
the gauge transformation. Now, regarding the Yukawa interaction, since only Nr., which 
is real symmetric matrix, appears in the 0{z,z) part of the wavefunction (4.82), all the 
theta function identities described in Sections 4.1, 4.2 hold for this new 0{z, z). Similarly, 
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as in the expression (4.49), the Yukawa couphng Yijk now has the following form, 



(|detN^,| 


IdetN^I) 




det N| 





E 



p,p 



1? 



detNi^detNl^ 



(0| (det det N^)^(N^"'n^N^- V) 

(4.95) 



with p running over all the states inside the cell generated by ei(ietN^N|j^ ^ and 
e2detN^'N^ ^ . Similarly p runs over all the states inside the cell generated by eldetN^N^ ^ 
and e^(ieiN^N^"\ 



4.8 Constraints on the results in section-4 and further general- 
ization 

To summarize, in this section we have given a close form expression for the Yukawa cou- 
plings in the magnetized brane constructions, when in general both oblique and diagonal 
fluxes are present along the branes. However, the results of this section are somewhat 
restrictive, since the basis wavefunctions used for the computations are well defined only 
when the intersection matrices satisfy a positivity condition given in eq. (2.17) for arbi- 
trary complex structure fl. A similar positivity criterion, for the case when p^ixj and Pyiyj 
are nonzero, can be written using the wavefunction (4.82), as well; it implies simply the 
positivity of Nr. 

On the other hand, in realistic string model building, one may need intersection ma- 
trices that are not necessarily positive definite. The simplest examples correspond simply 
to diagonal intersection matrices, having some positive and some negative elements along 
the diagonal. In such a factorizcd torus case, there is a unique prescription, to define 
the basis functions corresponding to the negative elements in the intersection matrix, as 
given in [5], consisting of taking complex conjugates of the wavefunctions for the positive 
elements. Such a prescription also works, in the case of oblique -|- diagonal fiuxes, when 
some intersection matrices are 'negative-definite' rather than being positive definite. One 
can then take a complete complex conjugation over all the coordinates, in order to obtain 
a well defined wavefunction. 

Such a process, however, does not work when oblique fiuxes are present and intersection 
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matrices have mixed eigenvalues. Note that a diagonal flux of the type F^i-i preserves its 
(1, l)-form structure, under the interchange : z'' — > z\ required by supersymmetry. This 
is, however, no longer true when oblique fluxes are present, since off diagonal elements of 
a (1, l)-form flux, say does not remain of the (1, 1) form when complex conjugation 

is taken only along or z'^. 

In order to cure the problem, one needs to construct new basis functions. We present 
the results of our investigation in the next section, where we first restrict to the case of a 

compactification, for simplicity. The complications arising from the oblique nature of 
the fluxes are manifest in the T'^ example as well, though it is possible to generalize the 
result to the full T^, which is discussed in Section 5.8. 

5 Negative-chirality fermion wavefunction 

As already mentioned, the basis wavefunctions given in cq. (4.3), used for deriving the 
Yukawa coupling expression in eq. (4.80), are constrained by the Riemann conditions 
(2.17), which imply in particular the positive-definiteness of the matrix N/mf2. 

Now, first restricting to T^, we will show that the basis function (4.3) corresponds 
to the positive chirality spinor on T^. On the other hand, to accommodate intersection 
matrices, having two eigenvalues of opposite signature, one needs to find out the basis 
function corresponding to negative chirality spinor. The need to use such basis functions, 
for intersection matrices with mixed eigenvalues, can be easily seen in the case when the 

factorizes into x and one turns on only non-oblique (diagonal) fluxes. In this 
case, the intersection matrix has one positive diagonal element along the first and one 
negative diagonal element along the second one. Good basis functions are then products 
of two wavefunctions of opposite chiralities [5] , and the total wavefunction on is of 
negative chirality. 

Our task therefore amounts to searching for the basis functions corresponding to neg- 
ative chirality spinors on with oblique fluxes. Search for fermion wavefunctions in the 
presence of arbitrary fluxes (in general oblique) has been pursued in [28]. However, the 
resulting wavefunctions are presented in terms of diagonalized coordinates and eigenvalues 

of fluxes. Any such solution is however unsuitable for the Yukawa computation, both for 
the purpose of extracting the selection rules of the type given in eq. (4.40), as well as 
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in actual evaluation, since the diagonalized coordinates become 'stack dependent' and in- 
herent nonlinearities involved in the diagonalization process appear in the wavef unctions, 
prohibiting the derivation of Yukawa couplings in a concrete form. 

In this section, we are able to write both the positive and negative chirality basis 
functions in a 'unified' fashion, by showing that all basis functions have a form similar to 
the one given in eq. (4.3). However, the complex structure Q appearing in cq. (4.3) for a 
positive chirality wavefunction needs to be replaced by an 'effective' modular parameter 
matrix fl = flU,, in order to accommodate the negative chirality wavefunctions, where 
is given in terms of the elements of the intersection matrices (as explicitly obtained later). 
We also show that our results reduce to the ones in [5] for the case of diagonal fluxes. 

First, in the next subsection we present new basis functions, relevant for the situation 
when the intersection matrices are neither positive nor negative definite. In a later subsec- 
tion, we show how the negative chirality spinor basis functions can be identified with the 
positive chirality ones given in eq. (4.3), with an effective modular parameter, defined in 
terms of the fluxes. We verify this fact by mapping the wavefunctions into each other, as 
well as, by showing explicitly that the relevant field equations transform into each other 
through such a mapping. As a result, we are able to absorb the complications associated in 
the diagonalization process of the modular parameter matrix, and the final wavefunction 
once again has an identical form as given in eq. (4.3), however, with a flux dependent 
modular parameter argument. 

5.1 Construction of the wavefunction 

In this subsection, as mentioned earlier, we discuss the case of 4-tori, though generaliza- 
tion can be analyzed in a similar manner. We flrst also restrict ourselves to the situation 
with canonical complex structure: Jl = i/2 and Q = i/3 for and respectively, where 
1(1 represents the (/-dimensional identity matrix. The generalization to arbitrary Q is given 
in subsections 5.6 - 5.8. Now, in oder to avoid the restriction to the positivity condition 
(2.17), we present an explicit solution of a wavefunction of negative chirality satisfying 
both the eqiiations of motion, as well as the periodicity requirements on T^. 

Going back to the positive chirality wavefunctions, note that the two equations for 
the component x+ ©Q- (^-H) (derived from the original Dirac equation (A. 6)) can be 
simultaneously solved, since when acting on x+ with two covariant derivatives, we have: 
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[Di, D2] ~ Fj^l and the RHS is zero, since all the (0, 2) components of the gauge fluxes are 
zero in order to maintain supersymmetry. The superscript ab in this relation implies that 
we need to take the difference of fluxes in brane stacks a and b due to the combination 
A" - that appears in eq. (A. 11) for the bifundamental wavefunction. Same is true for 
the two x+ equations, since (2, 0) components of the fluxes are zero as well. On the other 
hand, the relevant equations for the negative chirality spinors are: 

D^X- + D2X-^0, (5.1) 

and 

D2X- - D^x- = 0. (5.2) 

12 12 — 1 

When only one of the two components X- is excited at a time, X- satisfy: Dix_ = 
D2X- = or Dix^ = D2X- — 0. But none of these sets of equations can be consistently 
solved when oblique fluxes are present, since [Di, D2] ~ F12 ^ 0. 

1 2 

The two negative chirality components X- therefore need to be mixed up in order to 
obtain a solution of the relevant Dirac equations, when oblique fluxes are present. In other 

1 2 

words, we need to simultaneously excite both X- ■ Then, taking 

X- = aV', X- = P'^, (5.3) 

equations (5.1) and (5.2) become: 

{(3D2 - aD,)ij = 0, (5.4) 

and 

{pD, + aD2)iJ = 0. (5.5) 
In order for these two equations to have simultaneous solution, one obtains the condition: 

- a(3F^^ - a^F^l + + ajSF^^ = 0, (5.6) 

where F'^ = N,-^ is again the difference of fluxes in brane stacks a and b and N,-:; is 
the same hermitian intersection matrix, eq. (4.84), used in writing the positive chirality 
wavefunction and Yukawa couplings in eq. (4.3), and other parts of Section 4. When 
Px^x^ — 0, and ft — ils, N reduces to the real symmetric matrix as in eq. (2.15). 
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Fortunately, equation (5.6) has arbitrary solutions of the type: 

1 -q 



Tiab 



N = N- 



11 



q 



22 



-q q 



q 1 



(5.7) 



with q = ^ and Nii, N22 being arbitrary integers whose notation will become clear later 
(see eq. (5.26) below). The RHS of the above relation is a general parameterization of a 
2x2 symmetric matrix, since the two terms can be written as 



nab 



N = TV- 



11 





(5.8) 



After having shown the possible existence of the solution of the type (5.3), we pro- 
ceed to find the explicit form of the wavefunction ip by applying the allowed orthogonal 
transformations on the wavefunction of the negative chirality fermion on a which is 
factorized into x T^. To obtain the explicit form of this orthogonal transformation, we 
start by writing the coordinate coordinate, X'^ — z\z^ {i — 1,2), in the spinor basis. 
We note, for the choice of Dirac Gamma matrices (in a real basis) given in eqs. (A.l), 
(A.2) that 



M 



Zi Zi 



\ 



Z2 



Z2 
-Zl 



Z2 



-Zl 



(5.9) 



with Zi — Xi + ii/i and Zi — Xi — ii/i, {i — 1,2), which factorizes into 2x2 blocks providing 
the basis on which ,SC/(2)'s in the Lorentz group : SU{2)l x SU{2)r ~ SO{l, 3) act. We 
get in the spinor basis in the form of a 2 x 2 matrix: 




(5.10) 



Now to understand the transformation properties of the fermions on T^, we consider 
the following transformations on Xaa'- 




-j(6»i+6»2) 



Z2 



(5.11) 



Zl 
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We learn from eq. (5.11) that when factorizes into x T^, the transformations of 
the positive and negative chirahty fermions on the two T^'s can be read off from the 
transformation rules of zi and Z2 given above. Indeed, the transformation rules for the 
fermions ijj^^ on the two T^'s, denoted by indices i — 1,2 are: 

In this described in the section 3.1, the fermion wavefunctions can be written 

as a direct product of the ones on two T^'s as in eq. (3.3). We obtain the transformation 
of wavefunctions (eq. (3.3)): 

*i — > e^^^^i, — > e-^^2*L (5.13) 

It follows that a left transformation (^i 0,02 = 0) acts independently on (left handed) 
positive chirality wavefunctions, and a right transformation (^i — 0,02 ^ 0) acts on 
the negative-chirality (right handed) wavefunctions. Now, consider the following complex 
transformation on vectors in spinor basis: 



(5.14) 



Case-I: For e — h — 1, f — g — 0, c — —b, a — d, i.e a left transformation results in the 
following orthogonal coordinate transformation, 

Zl — > azi + hz2] Z2 — > az2 — bzi. (5.15) 

Case-II: Similarly, for a = d = 1, c = b = 0, h = e, f = —g, i.e a right transformation 
leads to 

zi — > ezi - fz2] Z2 — > ez2 + fzi. (5.16) 

In order to maintain the holomorphicity of the gauge fluxes, one therefore needs to 
make use of the later transformation, in order to generate a general wavefunction, starting 
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with the one which corresponds to the diagonal (non-obhque) flux. In addition, we need 
to maintain the integrahty of the fluxes, as we make such orthogonal transformations. 
However, in our case, we do not make use of any speciflc form of the transformation and 
rather use the above analysis as a guide for writing down a general solution. We then 
verify the equations of motion directly, in order to conflrm that the solution we propose 
is indeed the correct one. 



5.2 New wavefunction 

We now use the transformation (5.16) to obtain the wavefunction associated with the 
negative chirality fermion bifundamentals, starting with a wavefunction associated with a 
negative chirality spinor for a diagonal flux. In the notations of eq. (A. 5), it corresponds 
to exciting only the negative chirality component 

(5.17) 

We ignore the explicit form of '0, except to note that after the transformation (5.16), one 
generates 





(5.18) 



while (^+, ^+) remain zero. In the gauge sector, such wavef unctions are parameterized 
in the bifundamental representations by: 

XaA ^^^g^ 

as also given in eq. (A. 9). For negative chirahty components, the equations to be satisfled 
by the various components are: (see eq. (A. 10)) 

dix- + 52xi + {A' - A'),,x- + {A' - A'),,x- = 0, 

92X'- - dix'- + {A' - A%,x^_ - {A' - A%,x'- = 0. (5.20) 

We now show that the solution to eqs. (5.20), together with proper periodicity require- 
ments on T^, is given by the basis elements: 

^'^^^^ ■ f{z,z)-e{z,z) (5.21) 
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where, 

y^^^ _ ^in[{S!ijZiImzj)-{'NQZiImzj)] ^ (5.22) 

0(^^^^^^ ^ g7ri(i)[(mi+ji)M-(m3+jj)]g27ri[(mi+ji)N-2jg27ri^ (5.23) 

with 

M,j = N,j-Nij (5.24) 

where both N, N are real, symmetric matrices, given earher in eq. (5.7), and so also is 
M (Mjj — Mj-). We retain, however, both types of indices: i and j to incorporate real 
as well as complex components of the (1, l)-form fluxes Fjj. Also, an extra factor of i in 
the exponent of 0{z, z) corresponds to the fact that we are working with the canonical 
complex structure : D, — il2 ior the present example of the fermion wavefuncton on T^. 

The wavefunction (5.21) satisfies the Dirac equations (5.20) for the following gauge 
potentials: 

{A' - A%, = (Nn + Nii)zi + (Ni2 + Ni2)^2 

{A' - A'),-, = (Ni2 + Ni2)^i + (N22 + N22)^2. (5.25) 

The intersection matrix N is therefore given by: 

N = N + N, (5.26) 

as appearing previously in eqs. (5.6), (5.7). Also, we have imposed the following con- 
straints, in order to retain the holomorphicity of gauge potentials: 

a ^ -Nil ^ -N12 ^ N12 ^ N22 ^ 1 27) 
P N12 N22 Nn N12 1' 

Note that the ratios of the matrix elements of N and N arc identical to those given in eq. 

(5.7). We have therefore explicitly shown that the solution given in eqs. (5.21) - (5.23) 

satisfies the equations of motion. The transformation properties of this wavefunction (5.21) 

along the four 1-cycles of T^, are given by: 

^J,N,Nj-£_^ = gi7r([N-n]-Imi--[N-n]-Im2) . ^7,N,Nj-^^ 

(5.28) 

^J,N,N^-_^ ^ g-j7r([N-n]-ile2+[N-nl-Rel) . ^j,N,N^^^ 

provided that 
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• Nij = (N + N)fj e Z, i.e (N + N) is integrally quantized, 

• J satisfies: J ■ (N + N) e Z". 

We therefore notice that the integer quantization is imposed only on the intersection 
matrix N given in eq. (5.26) and does not necessarily hold for the matrix M in eq. (5.24). 
Explicitly, we have: 

N = N + N = Nii I ^ ~^ 1 + N22 I ^ ^ 

\-q q / \^ ^ 

where the first eq. in (5.29) is identical to the solutions in eq. (5.7). 

Note that the wavefunction given in eqs. (5.21), (5.22) and (5.23) is now well defined, 
as the series expansion in eq. (5.23) is now convergent. To show this, we note the following 
relation: 

detN = -detM = NnN22(l + 0^- (5-30) 

As a result, in the case when det N is negative ( when N has two eigenvalues of opposite 
signatures), detM > 0. So, the series (5.23) is now convergent when the two eigenvalues 
are of positive signature, since it is the quadratic part, in the summation index in theta 
series, that dominates in the exponent of this expansion. An overall complex conjugation 
will be required, for the case when two eigenvalues are negative rather than positive. 

5.3 Normalization 

Now that we have found a basis of wavef unctions, classified by the index ji in the exponent 
in (5.23), we proceed to show its orthonormality. The wavefunctions described in eqs. 
(5.21), (5.22), (5.23) can be re-written in terms of the real coordinates x and y as follows: 

^J,N,M _ _ gi7r[^-N-ir+iy-M-j/| ^ g7ri(i)[(m+j)-M-(m+J)]g27ri[(m+j)-N-£+i(m+J)-M-ji] (5 31) 

Then the following orthonormahty conditions are satisfied: 

f (V^^'N'M)*^J,N,M _ ^5 32) 
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To verify the orthogonality relation and obtain the normalization factor, we note that, in 
terms of the wavefunctions (5.31) we have: 



^^fc,N,M-j*^J,N,M _ j^^ . ^-in[x-TS!-y-iy-M-y\ ^■Ki{i)[{l+kyMil+k)] _ ^-2ni[{l+kyTSI-x-i{l+kyM-y\ 

feZ" 

J\f-, . ^inlx-TSI-y+iy-m-yl g7ri(i)[(m+J)-M-(m+j')] _ g27rj[(m+J)-N-x+i(m+J)-M-a] 

= A/W- • e"^'^^'^'^'^ g7ri(i)[(m+i')-M-(m+j')] _ g7ri(i)[(r+fc)-M-(r+fe)] 

g27ri[(m+j')-(r+fe)]-N-i _ g27ri(i)[(m+J)+(r+fc)]-M-jr gg^ 

The integration over x in eq. (5.32) imposes the condition j = k and equality on the sum- 

— * — * 

mation indices m — l. In particular, the condition j = k gives our orthogonality condition 
(5.32). One can now obtain the normalization factor by performing the integration: 

•1 

^-2TTy-M-y ^ g-27r(m+J)-M-(m+j) . g-47r(m+/)-M-j? 



meZ" 



/I 
'^^(^) g-27r((m+j)+^)-M.((m+j')+y) 
.meZ" 



(5.34) 



One can integrate over y, using 



LmeZ" J meZ"'^° 

oo 



J — t 



g-27r[(m+J)+j/]-M-((m+J)+j/)j 

(5.35) 



-27rj/'-M-2/' 



The integration (5.35) fixes then the normalization constant to 

= (2|detM|)'/' • Vol(T^)-^/2^ Vj. 

5.4 Eigenfunctions of the Laplace equation 



(5.36) 



The wavefunctions (5.21) not only represent zero modes of the Dirac operator, but are 
also eigenfunctions of the Laplacian. In order to see this, we start with computing the 
Dirac operator in four dimensions. In our notations: 

di 82 
82 -Bi 
V 82 -8, ) 



(5.37) 



49 



which leads to 



/ L>i A + D2D2 



\ 



DiDi + D2D2 



D2D2 + D,D, 



\ 



22 



DiDi + D2D2 J 



= A + 



—Fii + F22 



(5.38) 



22 



v 



(Fii + F22 



The Dirac equation — 0, with given in eq. (3.3), imphes that such basis functions 
are also eigenfunctions of the Laplacian A. The question whether massless scalars exist, 
depends on whether some combination of fluxes appearing in eq. (5.38) vanish. Of course, 
their existence is guaranteed in the supersymmetric case. 

5.5 Mapping of basis functions from positive to negative chiral- 



We now show that the basis for the negative chirality wavefunction, given in eqs. (5.21), 
(5.22), (5.23) can in fact be obtained by a mapping from the basis of the positive chirality 
wavefunction given in eq. (4.3). We also present the mapping between the corresponding 
field equations. Our mapping reduces to the ones in [5] for the case of factorized tori. 

More precisely, we show that our negative chirality wavefunction, given in eqs. (5.21), 
(5.22), (5.23), as well as (5.31) (for a trivial modular parameter matrix : Q = 2/2) is 
identical to the positive chirality wavefunction (4.3) for a 'nontrivial' (flux dependent) 
modular parameter matrix Q = iQ. Explicitly, Q is given in terms of the ratios (g) of flux 
components. This result gives a 'unified' picture of all the relevant basis functions. Later 
on, in Section 5.7, we show that a similar mapping holds for nontrivial complex structure 
on T^, by examining the equations of motion. 

Let us write down explicitly the wavefunction (4.3) for complex structure with arbitrary 



ity 
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fl (= iCl). 

^ ^ii7[x.lSI' .y+iCly.TSI' .y\ . ^ gi7r[(m+j)iN'n.(»n+J)]g2i7r[(m+J)(N's+iN'n.j?)] gg^ 

?neZ" 

where N is changed to N' to show a distinction between the two wavefunctions for the 
purpose of defining the mapping as given below. Next consider the negative chirahty 
wavefunction (5.31), written in terms of real coordinates x and y, 

^J,N,M ^ ^in[x-'N-y+iy-'M-y\ g7rj(j)[(m+J)-M-(m+J)]g27rj[(m+J)-N-J+i(m+J)-M-y] ^q-j 

mSZ" 

It is now easy to check that the above equations (5.39) and (5.40) precisely match with 
the following identification : 

N = N + N = N', 

M = N-N = N'0^0 = N-^M, (5.41) 

with Q = ifl, and is a real matrix. For the N and M, defined in eq. (5.29), and fl 
are given by; 



:i+g2)2 



1 -q\ I { q^ 1 



Nil \-q q^ N22 \ g 1 



(5.42) 



(l + ?2) ^ ^ ' 

We have therefore shown explicitly that the positive chirahty basis wavefunction (4.3), 
known earlier in the literature, can be mapped to the negative chirahty wavefunctions 
that we have constructed in eqs. (5.21)-(5.23), (5.31). Such a map also confirms the 
validity of our construction for the negative chirahty basis functions, presented using basic 

principles, such as equations of motion as well as periodicity requirement. In fact, in the 
next subsection, the same mapping is also obtained through comparison of the relevant 
equations of motion, which further confirms our results for the construction of the basis 
functions. Note that for g = or g — > 00, corresponding to the case when both matrices 
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N and M in eq. (5.29) are diagonal, we have: 

/ -1 n\ 

(5.44) 





(i o) 




(-1 o\ 






1 , or Q = 1 






[0 -Ij 




10 1 



respectively. As a result, one reproduces the known mapping of the wavefunctions between 
positive and negative chirality spinors in the case when is factorized into x [5] . 

5.6 Mapping the equations of motion 

In order to derive a similar mapping of the equations of motion, we show below that the 
covariant derivative operators appearing in eqs. (A. 12) for the positive chirality wavefunc- 
tion, with a nontrivial complex structure [ifl), are equivalent to the derivative operators 
appearing in eqs. (5.4), (5.5) for the negative chirality wavefunction (with complex struc- 
ture O = i/2). The mapping of corresponding gauge potentials can also be shown in the 
same manner, since they have similar dependence on the complex structure as the deriva- 
tive operator. Note that the complex structure appears in the wavefunctions as modular 
parameter matrices. We therefore reconfirm the mapping between the two wavefunctions 
by comparing the equations of motion as well. 

We now examine the Dirac equations for both cases. For the first one, with arbitrary 
Q(= iCl), we have 

Z + Z _ ,Av_i ( Z — Z 



x + iCly; z — X — iCly =^ x — ""^ ; y—{d) 



2i 



which implies 



dzi ~ 2 \dx' dyi 

d If d .,^,1 d 



Then, the Dirac equation for the positive chirality wavefunction is: 

D,^^^^'''{z,n) = ^ (d,. + z(0)t;d,,) ^^''''{z,n) = 0, i,j = 1,2. (5.46) 

On the other hand, for the negative chirality solution (5.20), with complex structure 
fl = il2, the relevant derivative operators are: 

(/3Di + aD2) ^^'^'^ = 0; - aD^) ^^"'^'^ = 0. (5.47) 
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These equations, using the definitions z'^ — -\- iy\ Zi = — iy\ can be rewritten as: 

Now using ^ = q from cq. (5.27) and comparing the equations (5.46) and (5.48), one finds 
that they precisely match for the following complex structure: 

(1 + q') \-2q q^-lj 

which is exactly the same as eq. (5.43). Thus, the wavefunctions as well as the Dirac 
equations for both cases match exactly. This mapping can be generalized further, as given 
in subsection 5.8 below. 



5.7 Mapping for arbitrary complex structure Q 

In this subsection, we generalize the mapping between the equations of motion associated 
with the positive and negative chirality wavefunction to the case of compactification 
with arbitrary complex structure Q. Now, the negative chirality basis functions satisfy: 

^ + >(S1)- (^jfl^B,,) + .(fi)- (^^C,.) } '/^'•''•" = , (5.50) 

which can be identified with the equations satisfied by the positive chirality wavefunction 
with O = 00, as can be seen through the decomposition: 

dzi 2\dx^ dyi) ' 

Thus, eq. (4.3) with Q = QQ, with Q given in eq. (5.49), provides the negative chirality 
solution for arbitrary complex structure O, where both 'oblique' and diagonal fiuxes are 
turned on. 
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5.8 Generalization for the T^- case 

In this subsection, we generalize the results obtained so far for negative chirality fermions 

on T'* to the more general case. Wc only consider the wavefunctions that are well 
defined with two positive and one negative eigenvalues of the 3x3 Hcrmitian intersection 
matrices, since these will complete the list of well defined wavefunctions, once complex 
conjugations are taken into account. For the case of T^, the relevant equations, obtained 
by generalization of eqs. (5.4) and (5.5) to be examined, are: 

{aDi - (3iDi)^ = 0, (5.52) 

and 

{aDi + (3iD^)iP = 0. (5.53) 

Note that in these equations and below, the indices i,j — 1,2 (used for the with wave- 
functions of positive chirality) . In order for the above two equations to have simultaneous 
solution, one obtains the condition : 

a^i^f + aPiF^^ - a(3jF;f - MjF^f = 0, (5.54) 

where F""^ = N is the difference of fluxes in brane stacks a and b. The general solution of 
this equation is of the following type: 

= N = AT vy; ^ vy; y y ^^^^ 

\-q q{q)^J \ Nq N ) 

where TV is a 2 x 2 matrix and TV is a number. Also, q is the two-dimensional {2d) vector 
defined as: 



(5.56) 



with qi = ^. 

Now, after showing the possible existence of the solution by defining F"'' in (5.55), for 
the negative chirality wavefunction on T^, we proceed to present a mapping between the 
equations of motion for negative chirality and positive chirality wavefunctions on T^. As 
described before in section 5.6. Here also we show that the covariant derivative operators 
appearing in eqs. (A. 12), for the positive chirality wavefunction, with a nontrivial complex 
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structure are equivalent to the derivative operators appearing in eqs. (5.52), (5.53) for the 
negative chirahty wavefunction (with complex structure Q = i/3) and the corresponding 
gauge potentials map in the same manner. 

For the positive chirahty case, with arbitrary Q{— iQ) and eqs. (5.45), (5.45), the 
Dirac equation reads: 

D-^J^'^^'iz, n)^^ [d,, + mzlDy.) V''>'(z, Q) = 0, 1/ = 1, 2, 3 . (5.57) 

On the other hand, for the negative chirahty solution, with complex structure = i/3, 
the relevant derivative operators, given in eqs. (5.52), (5.53), take the form: 

^ { {a^Sij + A/9,) D,, - I (2aA) Dy. + i {^i^j - a^S^j) Dy, ] V^^^'^'^ = 0, 

^ { {Pi + a") +t{a^- (3f) - t {2(3,a) Dy. ] ^^'^ = 0. (5.58) 

Now, defining new 2x2 matrices, 
and 

Pi = (2aA), 

eqs. (5.58) can be re-written as: 

i {O., - i (A-'P), D,, + i (>l-'B)y D,, } = 

A comparison of equations (5.57) and (5.60) implies that they precisely match for the 
following complex structure: 

This expression for the complex structure generalizes the one derived earlier in eq. (5.43) 
for the case. The results are also easily generalizable to arbitrary complex structure Q 
following the discussions in subsection 5.7 for the special case of (see eq. (5.51)). 
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5.59 



5.9 Computation of Yukawa couplings 

Now that we have derived both the fermionic and bosonic internal wavefunctions and 

expressed them as an orthonormal basis, we compute the Yukawa couphngs using the 
basis wavefunctions (5.31). We also point out how the results derived below reduce to the 
ones in section 4. 

Starting with basis functions described in eq. (5.31), for the case of the canonical 
complex structure Q — ih (in the case), we have: 



^?,Ni,Mij---^ _ ^i,N2,M2,^-j ^ _ gi7r[x-(Ni+N2)-i?+jy-(Mi+M2)-j^ (5.62) 



^7ri(i) [(/i +?) -Ml • (h +?) + (^2 + J) • M2 • (^2 + J)] 



g27ri[(Zi+i)-Ni+(Z2+J)-N2]-ig27ri(0[(/i+?)-Mi+(/2+J)-M2]-j? 



This expression can be re- written as: 

^i,Ni,Mij-^-j _ ^i',N2M2,|-^j ^ _ _ gi7r[£-(Ni+N2)-r7+ij7-(Mi+M2)-{;i (5.63) 

■ e' 



3 

g7ri(i) (f-Q-l) g27ri(f -Q-X) _g27ri(i) (f -Q- Y) 



ii,/2ez" 



where we defined the 4d- vectors: 






, Y= M , (5.64) 



and the 4(i-matrices: 

/ Ni \ . / Ml , ^ ^ 

Q= 'Q= ■ 5.65 

\^ N2 y \^ M2 

Using the transformation matrix T, defined in eq. (4.9), and eqs. (4.10)-(4.14), we 
explicitly write the terms appearing in the exponents in the RHS of eq. (5.63) as: 

(1)^ • Q • (5 = {If . (T-'T) ■ Q • (T^iT-Y) ■ (f), 
(f • Q • X) = f • {T-^T) ■ Q • {T^{T-Y) ■ X, 
(f • Q • Y) = f • (T-^T) ■ Q • (T^(r-i)^) • Y. (5.66) 
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Then using: 

^ /(N1 + N2) , ^ ^ 



(Ml + M2) (MiNj-^ - M2N2~^)q;^ 

Q;(Nr^Mi - N2"^M2) Q;(Nr^MiNr^ + N2~^M2N2-^)q;^ 



.,T^-i , (i + M(Ni-^ + N,-^)-iN,-^ + (j+/2)(Nr^+N,-VNr^ / , _ 

(1) ^ = r.^ . , ' (5-68) 



and 



(a-T(Nr^ + N2-^)-M(?+/0-(J+y] I' ^-^ 



we can re- write eq. (5.63) as 

^?,Ni,Mi , ^l,^2M2,^g^ ^ . . gi7r[i'-(Ni+N2)-j?+if(Mi+M2)-y] ^ (5-71) 
g7ri(i)({[(;i+?)Ni+(i;+J)N2](Ni+N2)-i}-(Mi+M2)-{(Ni+N2)-i(Ni(i+Zi)+N2(J+i;^ ^ 

g27ri{[(Zi+?)-Ni+(i;+i)-N2](Ni+N2)-i}-(Ni+N2)i- . g27ri(i){[(iI+?)-Ni+(i2+J)-N2](Ni+N2)-i}-(Mi+M2)j/ ^ 
g27ri(0{[(i+«0-(J+i2)](N-i+N-i)-ia-i}-a(Ni-iMi-N2-iM2)-ir ^ 

g7ri(i)({[(«7+i)Ni+(i2+J)N2](Ni+N2)-i}-(MiNi-i-M2N2-i)a^{(a-i)^N2(Ni+N2)-iNi[(?-J)+ ^ 
g7ri(i){[((i-i)+(Zi-Z2))Ni(Ni+N2)-lN2a-l]-[a(Ni-lMi-N2-iM2)](Ni+N2)-MNi(i+/^ ^ 
g7ri(i){[((?-J)+(ii-Z;))Ni(Ni+N2)-iN2a-il[a(Ni-iMiNi-HN2-iM2N2-i)a^][(a-i)'^N2(Ni+N2)-iNi[(^^^^ 

Now, in a similar exercise as the one performed earher in sections 4.2, 4.3, 4.4, we 
rearrange the series in eq. (5.71) in terms of new summation variables I3, U, m, whose values 
and ranges are assigned as in these sections.^ With the value of a = (detNi detN2)/, 
^For details see sections 4.1, 4.2, 4.3, 4.4. 
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defined in eq. (4.21), eq. (5.71) takes the form: 

^?,Ni,Mi _ ^J,N2M2,j^^^ ^ j\J-_^ _ j\J-_^ , gi7r[i-{Ni+N2)-?+i?-(Mi+M2)-j?] (5-72) 



"S^ g7ri(i)[(iNi+jN2+mNi)(Ni+N2)-Hi3]-(Mi+M2)-[(Ni+N2)-i(Nii+N2i+Nim)+i3] ^ 



g27rj[(?Ni+jN2+mNi)(Ni+N2)-Hi3]-(Ni+N2)5 . g27ri(i)[(?Ni+jN2+mNi)(Ni+N2)-i+i3]-(Mi+M2)j? ^ 
g27ri « [(?- j+m) ^il|^^l±^2^^ 

^7ri(i)[(iNi+jN2+mNi)(Ni+N2)-i+i3]-[{detNidetN2)(MiNi-i-M2N2-i)]-[5^^^^^ ^ 

Using from eq.(5.31): 

^ ^ g7ri(i)[(z|+fc).M3-(«^+fc)] . g-27ri[y7+fc).N3-^-ia|+fc)-M^ ^g^^g-j 

we can then proceed to calculate the Yukawa coupling: 

Yijk = C7abc9 f dzdzi ■ V'^'^^'^^ • V''''^^^^ • (V^*^'^-^«)* {i = 1, 2) . (5.74) 
Consider first the integration over x: 

J -gj7r{^-[{Ni+N2)-N3]-j/} ^ ^ g27ri[(iNi+jN2+?nNi)(Ni+N2)-i+/3]-(Ni+N2)ig-27ri(;|+fc)-N3-x 

(5.75) 

which implies, using (Ni + N2) = N3 , the following conditions: 

• equality of the summation indices I3 = l'^, 

• the relation (ilMi + JN2 + mNi)(N3)-^ = k . 

Note that (Ni + N2) = N3 is a valid condition in a triple intersection since lab + Ibc — lac, 
with complex conjugation taking care of the fact that lac — —Ica, which changes the signs 
of N3 and M3. Also, as in section 4.3, 4.4, for any given solution of the above constraint 
equation for i, j, k, rh, other solutions inside the cell of eq. (4.34) that are shifted by m's 
satisfying mNiNa"^ : integer are also allowed. In view of this, as in eq. (4.47), we break 
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the sum over m into two parts, one corresponding to m, which is a given specific solution 

— * 

of eq. (4.40) and the other ones as given by sum over integer variables p and p whose 
ranges are as defined in eq. (4.45). 

Imposing the constraints from the x integration, we obtain: 

Yijk^Uabcg-MrM^-N-^ (5.76) 



j (^2-|g-7r[jr.(Mi+M2+M3)-S] ^ ^ 



g7ri(i)[fe+/3]-(Mi+M2)-[fc+i3] X 



g7ri{i)[fc+/;].[(detNidetN2)(MiNi-i-M2N2-i)]-[5^2^5l±^|^^ ^ 



TTl' 

e 



^ gTTi « [fe+Zg] • (Ml +M2) . g^ri (i) [(i- J+m) ""I'^j^+ll^Nl""^ +«4] • [(det Ni det N2) (Mi Ni "i -M2N2 " ^ )] | 

where the range of the sum over p,p is as used in eq. (4.45) in section 4.3. 

The above expression for the Yukawa interaction can be written as following: 



Yijk = (7„6e^-(2=^)i(|detMi|.|detM2|.|detM3|)i (yoZ(r^)) 

^2 -|g-7r[jr-(Mi+M2+M3)-i/] ^ ^ g7ri(i)[K+L]-Q'-[K+l 

(^afecC/ • (2^)^(|detMi|.|rfetM2|.|rfetM3|)2 {Vol{T^))~^^ 



j ^2 -|g-7r[j?.(Mi+M2+M3)-j/] . ^ 



where we defined new 4(i-vectors: 



K 







(Y'lzQ')} (5.77) 




A; 

ni J ^ rn)\[ HetNidetN-^ J 



(5.78) 



(Mi + M2)y 
[(detNi detN2)(MiNi-^ - MaNa"^)] • y 



(5.79) 



and the 4(i-matrix: 



(M1 + M2) (detNidetN2)(MiNr^ -M2N2-')] 

l(det Ni det N2) (Nr^Mi - N2~^M2) (det Ni det N2)2(Ni-^MiNi-^ + N2~^M2N2~^)J 
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(5.80) 



— * 

with k appearing in eq. (5.78) restricted by the Kronecker delta relation written above, 



as following from the x integration, in eq. (5.75) and the range of the sum over p,p is as 
used in eq. (4.45) in section 4.3, we skip the details regarding them. 

In fact, the form of the result (5.77) is valid for all basis functions, whether correspond- 
ing to positive or negative chirality wavefunctions, since the negative chirality wavefunction 
(5.31), written for the complex structure fl — i/2 and used in obtaining the final answer 
for Yukawa coupling in eq. (5.77), reduces to the one for positive chirality wavefunction 
for the same complex structure when M is set to N (see eq. (4.3) for the general form 
of the positive chirality wavefunction). For such a choice: Mi = Ni, Q has a factorized 

— * 

block form and the vector Y' in eq. (5.79) now has a form: 



The theta function in eq. (5.77) then factorizes and the final answer reduces to the form 
given in eqs. (4.39), (4.49) for the choice r — i corresponding to the complex structure of 
our choice in the negative chirality wavefunction (5.21). 

The Yukawa coupling expression (5.77) can be further generalized to other situations. 
First, although the above analysis was very specific to the case of T"^ due to our choice 
of wavefunction in eq. (5.31), the generlization to the is staightforward. Mapping 
between matrices N and M is identical and follows from the definition of 0, in subsection 
5.8. The final answer is identical to the one given in eq. (5.77). 

Further generalization to the situation of arbitrary complex structure should also be 
possible, using the wavefunctions that emerge due to the mappings obtained in subsection 
(5.7) and scaling procedure presented in section (4.6) for the positive chirality wavefunc- 
tions. One, however, also needs to examine the symmetry property of the matrices Ni7r2 
etc., appearing in the definition of the wavefunction. We leave further details for future 
work. 

6 Mass generation for non-chiral fermions 

In this section, we briefiy discuss one of the applications of the results derived in the paper, 
for giving mass to the non-chiral gauge non-singlet states of the magnetized brane model 




(5.81) 
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discussed in [27]. This is a three generation SU{5) supersymmetric grand unified (GUT) 
model in simple toroidal compactifications of type I string theory with magnetized D9 
branes. The final gauge group is just SU (5) and the chiral gauge non-singlet spectrum 
consists of three families with the quantum numbers of quarks and leptons, transforming in 
the 10 + 5 representations of SU{5). Brane stacks with oblique fiuxes played a central role 
in this construction, in order to stabilize all close string moduli, in a manner restricting 
the chiral matter content to precisely that of SU{5) GUT. Another interesting feature 
of this model is that it is free from any chiral exotics that often appear in such brane 
constructions. However, the model contains extra non-chiral matter that is expected to 
become massive at a high scale, close to that of SU (5) breaking. 

The results of the previous sections can be used for examining the issue of the mass 
generation for these non-chiral multiplets in a supersymmetric ground state. The aim is 
to analyze the D and F term conditions, and show that a ground state allowing masses 
for the above matter multiplets is possible. The exercise will further fine tune our SU (5) 
GUT model to the ones used in conventional grand unification. 

Although, we will not be evaluating any of the Yukawa couplings explicitly, which using 
our results of the previous sections is in principle possible to do, the aim of the exercise 
below is to show that indeed one can give masses to non-chiral matter. Our procedure 
involves the analysis of both the F and D-term supersymmetry conditions. In the context 
of our previous work [27] , we like to remind the reader that certain charged scalar vacuum 
expectation values (VEVs) were turned on in order to restore supersymmetry in some of 
the "hidden" branes sector. These charged scalar VEVs gave a nontrivial solution to the 
D-term conditions, but left the F-terms identically zero in the vacuum. In the following, 
on the other hand, our aim is to find out the possibility for a large number of scalars in 
various chiral multiples to acquire expectation values. For this, we need to examine both 
the F and D conditions, as already mentioned. 

The model in [27] is described by twelve stacks of branes, namely C/5, C/i, Oi . . . , Og, A, 
and B. The magnetic fiuxes are chosen to generate the required spectrum, to stabilize all 
the geometric moduli and to satisfy the RR-tadpole conditions as well. The fiuxes for all 
the stacks are summarized in Appendix C. The fiuxes for stacks U^jUi, A, B are purely 
diagonal whereas stacks Oi . . . , Og carry in general both oblique and diagonal fluxes. All 
36 closed string moduli are fixed in a A/" = 1 supersymmetric vacuum, apart from the 
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dilaton, in a way that the T^-torus metric becomes diagonal with the six internal radii 
given in terms of the integrally quantized magnetic fluxes. 

The two brane stacks C/5 and Ui give the particle spectrum of SU{5) GUT. We solve 
the condition lu^Ui + lu^u^ = — 3 for the presence of three generations of chiral fermions 
transforming in 5 of SU (5) and continue with the solution lu^Ui = 0, Iur,Uj = ~3. The 
intersection of C/5 with Ui is non-chiral since lusUi vanishes. The corresponding non-chiral 
massless spectrum consists of four pairs of 5 + 5, which we would like to give mass^. 
Obviously, we would like to keep massless at least one pair of clcctroweak higgses but 
this requires a detailed phenomenological analysis that goes beyond the scope of this 
paper. Here, we would like only to show how to use the results of the previous sections in 
order to give masses to unwanted non chiral states that often appear in intersecting brane 
constructions. 

So, we have the following non-chiral fields where the superscript refers to the two 
stacks between which the open string is stretched and the subscript denotes the charges 
under the respective [/(l)'s :(0+'L^\0^'!j^\ 4), with numbers in the brackets denoting the 
corresponding multiplicities. Similarly, the intersections of the U5 stack with the two extra 
branes A, B and their images are non-chiral, giving rise to the extra 5-1-5 pairs: 
149), {(t)^+f\(t)^jf, 146), {(t)^'^ ,(t)^J^ , 51), {(t)^+f\(t)^J^\ 16). A common feature of all these 
states is that they arise in non-chiral intersections, where the two brane stacks involved 
have diagonal fluxes and are parallel in one of the three tori. It is then straightforward 
to give masses by moving, say, the C/5 stack away from the others along these tori. In the 
language of D% branes, this amounts to turn on corresponding open string Wilson lines. 

On the other hand, analysis of the particle spectrum on the intersections of the stack 
f/5 with the oblique branes Oa and O* , satisfying the condition lu^a + lu^^a* = 0, for a = 
1, .., 8 , leads to 4 x (23 + 14) = 148 pairs of (5 + 5) representations of SU{5): 

Iu,o. = -23, /[/,o; = 23, a=l,...,4, (6.1) 

Iu,o. = -U, /c/50s = 14, a = 5,...,8. (6.2) 

We then have the following chiral multiplets, 23), {(!)++% 23), ((/)?+^^ 14), (0++°% 

14) (a = 1, . . . , 4, 6 = 5, . . . , 8). In order to examine the mass generation for these fields, 
one needs to write down the superpotential terms involving the above chiral multiplets, as 
''For details, see Section 3.7 of [27]. 
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well as those coming from the brane stacks Oi, • • • ,08 and their orientifold images. The 
list of the later, involving purely oblique stacks, is given in Appendix C. 

Now, using the results in Appendix C in eqs. (C.13) and (C.14), one can analyze the 
associated superpotential and D-terms and look for supersymmetric minima. The relevant 
superpotential reads: 

= E (^i-^y ^^+?y ^^-^-'f + E ^^++''y 

ijk ijk 

+ E ^^i-'y (^"if'^y + E (^'^-'y (^+f'y ^^--^f (^-s) 

ijk ijk 

+ E ^ff {€'-^'y {4>if'y {4>-''-'f + E ^^i-^y ^^+?y ^^-^ 

ijk ijk 

where the sum over i,j,k runs over the "flavor" indices. The couplings Wq.^, given in eq. 
(6.3), can be read off from our results in the previous sections. In addition to the complex 
structure, these also depend on the first Chern numbers of the branes in each triangle. 

The F-flatness conditions (Fj) = {V^^W) — (at zero superpotential, = 0), imply 
that for each "triangle" at least two fields must have a zero VEV in order to form a 
supersymmetric vacuum [30]. In this theory, there exists indeed a supersymmetric vacuum 
where six charged fields remain unconstrained by the F-flatness conditions. Let's choose 
them to be {(fj^'), {<fj^'), (0-i°'), (0-!°'), (0^i^'), {(p-') (they are neutral under 
the U (1) of the U (5)). The remaining fields appearing in the superpotential acquire a mass 
from the F-term potential only if these unconstrained scalars possess a non- vanishing VEV. 
Indeed, their masses read: 



Ml 


-Ml , 






Ml , 


~ Ml 






Ml 


~ Ml 







(l'^o^o'/),M|_~M|J^~(|0„.„J^), (6.4) 

where ^nsof, denotes linear combinations of 0„go7 with 0^503 and <^oio'^ denotes linear com- 
binations of 00^07 with 00^03 ■ Thus, the leftover massless states from the intersection of 
C/5 with the oblique branes are 60 pairs of 5 -|- 5: 0U5O* for a = 1, 2, 5 of positive chi- 
rality together with the negative chirality states 0u6Oa for a = 6, 7, as well as 23 linear 
combinations of 0^503 with 0^507, and 14 0^504 ■ 

However, switching on non-zero VEVs for these fields, modifies the existing D-term 
conditions for the stacks of branes Oi, ....Og. Recall that, in [27], the stacks C/5, 0\. . .0^ 
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satisfy the supersymmetry conditions in the absence of charged scalar VEVs, but VEVs 
for the fields (j)^^, 4'+^+* 0+- switched on, for the same supersymmetry to be 
preserved by the stacks C/i, A and BJ The D-terms for each C/(l) factor of the eight 
branes Oi, Og read 



Do, = -|(/)^i^3|2 _ |0O3O*|2^ ^ -\(j)020l\2 _ |0O4O?|2 

Do, = -\(l)0*0^2 _ |0O.O||2 ^ _ -|0O3OS|2 _ |0O.O||2 



(6.5) 



We can regain the supersymmetry conditions Da = 0, Va = 1, . . . , 8 with C,a{F"', J) = 0, 
by switching on VEVs for the following fields: {(f>+f'), {(p+f'), {(p+f'), {<P+f'), {<P+fb, 
(0++^®), provided these fields do not modify the superpotential (6.3). The modified D- 
terms read: 



Do, = 

Do, = 

Do, = 

Do, = 

Don = 

Doe = 

Dor = 

Dos = 



^0^0^2^ 
^020^2^ 
0OiO3|2_ 
0O2O^|2_ 

0O5Oe*|2 + 

0O4O?|2 _ 
0O3O||2 _ 



0OiO,*|2 
|2 

0O4O?|2^|^O3Oa|2^ 1^040,* |2 
0O6O*|2 

0OvO*|2^|^O.O;|2^|^O3O*|2 
^OtOJ 1 2 ^ I^OeOg* |2 ^ |2 



(6.6) 



The supersymmetry conditions Da = 0, Va = 1, . . . , 8 with ia{F"', J) = can be simulta- 
neously satisfied if and only if the VEVs for all these fields appearing in the expressions 
(6.6), have the same value, say v^. Moreover we can restrict v « 1 in string units, as 
required by the validity of the approximation for inclusion of charged scalar fields in the 
D-term. 

We have therefore shown the mass generation for a large set of non-chiral fields as given 
in eq. (6.4). It is possible, that remaining ones can also be made massive by incorporating 
non perturbative instanton contributions to the superpotential. However, we leave this 



^For details see Section 5 of [27]. 
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exercise for the moment. We also do not give any superpotential couplings, in terms of 
fluxes, as given explicitly in the previous sections. 



7 Discussions and Conclusions 

In this concluding section, we first comment on the case of magnetized branes with higher 
winding numbers. The form of the wrapping matrices [13] for D9 branes on was 
discussed in our earlier papers [12,27]. They are real 6x6 matrices giving the embedding 
of the brane along spatial internal directions. The situation where worldvolume coordinates 
are identified with the spatial coordinates corresponds to W being diagonal. Then, for 
example, for a canonical complex structure — ils, the spatial components of the flux 
matrices are of the form given in eqs. (B.3), (B.4), (B.5). Taking into account the gauge 
indices, one obtains a block diagonal matrix structure for the fluxes, that reduces in the 
case of factorized tori to the form: 



winding numbers along various 1-cycles of {TY e T^- Also, N" and N'' are the number 
of branes in stacks a and b respectively and the above expression has a straightforward 
generalization when many such brane stacks are involved. 

In [5] , a gauge theoretic picture of the magnetic fluxes along brane stacks with higher 
winding numbers (> 1) was given. For instance, consider the simplest choice A'^" = A'^'' = 
1. In this case, the configuration of the brane stacks a and b with one Z^-brane each, 
having wrapping numbers n"', and 1st Chern numbers m"-, m^, is given by a flux matrix 
associated with a U (n"+n'') gauge group with flux having the internal (gauge) components: 



along the i'th and m^,n^ etc. are relatively prime. 

Given the C/(n" + n^) flux in eq. (7.2), the fermion wavefunctions associated with 
bifundamentals were constructed in [5]. The new feature is that, to have proper periodicity 




(7.1) 



with a and b representing the brane-stacks and i denotes the i'th T^. Also m'^''^ are the 
first Chern numbers, as given in eqs. (B.3) and (B.4), whereas n"''' are the product of the 




(7.2) 
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property for these fermion wavef unctions, non-abelian Wilson lines need to be turned 
on. In turn, these non-abelian Wilson lines mix up nf x components and the set of 
periodicity constraints only allows the bifundamentals belonging to the representations of 
the gauge group: U{P^) x U{P^), with = g.c.d.{m1, n"). In our example above we have 

pa _ p6 _ 1 

i i 

The case of obhque fluxes brings in extra complexities in the analysis due to the 
presence of six independent 1-cyclcs along which non-abelian Wilson line actions need 
to be fixed. Given the action of these Wilson lines, one can then proceed to obtain the 
wavefunctions as well as the Yukawa couplings. However, unlike the factorized situation in 
[5] , one finds that the action of non-abelian Wilson lines on the wavefunction, is dependent 
on the particular model, or more precisely, on the details of the oblique fluxes that are 
turned on. Further analysis along this line is, though cumbersome, possible. We now 
conclude our paper with the following remarks. 

In this work, we have been able to explicitly generalize the Yukawa coupling expressions 
to the situation when the worldvolume fluxes that arc responsible for moduli stabilization, 
chiral mass generation, supersymmetry breaking to iV = 1 etc., do not respect the factor- 
ization of into (T^)^. For the factorized tori, the mappings of the Yukawa couplings, 
superpotentials and Kahler potential between the type IIB and IIA expressions was dis- 
cussed in [5]. In the IIA case, the results are obtained through a 'diagonal' wrapping of 
the D6 branes in three T^'s. 

It will be interesting to map our IIB expressions, given in this paper to the IIA side 
and find the corresponding intersecting brane picture. As stated earlier, such a IIA con- 
struction will require putting the branes along general SU{3) rotation angles and then 
obtain the area of the triangles corresponding to the intersections of three branes giving 
chiral multiplets. 

Supersymmetry breaking is of course an important issue in model building. Though 
generally, for magnetized branes, one encounters instabilities in such a situation, it should 
be however possible to obtain non-supersymmetric magnetized brane constructions for a 
rich variety of fluxes accompanied by orientifold planes which can possibly project out 
tachyons that may be generated during the process of supersymmetry breaking. 

The recent developments in writing the instanton induced superpotential terms are 
also encouraging, for the purpose of examining the supersymmetry breaking as well as 
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up-quark mass generations in a GUT setting. In this context, it has been shown that 
the magnetized branes too can give rise to interesting superpotentials through the hft of 
fermion zero modes when fluxes are turned on. 

Recently, there have been interesting developments in deriving particle models and 
interactions from string theory, using the non-perturbative picture of F-theory [31-38], 
with a geometric picture of 7-brane intersection curves inside del Pezzo surfaces giving the 
chiral spectrum as well as Yukawa interactions, including those of the spinors of 5*0(10) 
GUT, and thus generating observable masses for both up and down type quarks. It is 
also interesting to note that F-theory results are reproducible in a globally consistent IIB 
string theory, taking into account the instanton generated superpotential terms [39]. It 
will be of interest to see the implications of these results on the construction of SU{5) 
GUT in [27], as well as on the Yukawa interactions discussed in this paper. 

Finally, it will be interesting to explore the generalization of our results to higher-point 
functions (computing couplings of higher dimensional effective operators) [40] and make 
explicit comparisons of our results with those in [14,15], where the situation with diagonal 
intersection matrices Ni, but non-factorized complex structure, is addressed through a 
computation of twist field correlations. However, one then needs to examine the effect of 
supersymmetry conditions (2.7) and (2.8) to see if the interaction indeed remains nontrivial 
in a supersymmetric set up. 

We hope to return to all issues above at a later stage. 
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A Wavefunction 



Wc first present the construction of chiral fermion wavefunctions on tori and give their 
representation in terms of theta functions. For definiteness we first discuss the case of 
4-tori, though chiral multiplet structure can be analyzed in a similar manner. To be 
explicit, for the moment we restrict ourselves to the canonical complex structure: Q = i/2 
and Q = i/3 for and respectively, where la represents a d-dimensional identity 
matrix. The general complex structure is restored while writing the wavefunctions as well 
as interaction vertices. 

To obtain the Dirac wavefunctions in T^, we start by writing four Dirac Gamma 
matrices (in a complex basis) : 

/O 2 \ / 2 0\ 





-2 

V oy 



^ / X cr" 



2 



\o y 



(A.l) 



where the information about the complex structure in the above expression is hidden in 
the fact that we have used the definitions: Zi — Xi + iyi in writing these Dirac matrices. 
Similarly, 

/o \ / o\ 



r^"i = a' X 



2 




V -2 0/ 
They satisfy the anti-commutation relations: 



, = / X a' 





2 
\0 2 J 



(A.2) 



|r^i , r^^ } = 0, {T'^ , T'^ } = 0, {r^^ , r^^ } = 45,, 

with i,j = 1,2. In the above basis takes the form: 

/l \ 
-1 



V 



(A.3) 



(A.4) 
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with 4-component Dirac wavefunctions having the form: 

^2 



(A.5) 



In such a decomposition of ip, Dirac equations for fermions in the adjoint representation 
are of the form: 



0, 
0, 
0. 



8,^1 + 82^1 + [A,„^^] + [A,,, 
82^1 - d,^[_ + [A,,,^V^ - [A,,,^\ 



(A.6) 



In a generic model, chiral fermions arise either from the string starting at a brane 
stack-a and ending at another brane stack-6 (or its image h*) or from strings starting 
at a brane stack a and ending at its image a* . We aheady showed the correspondence 
between a stack of magnetized branes and flux quanta in supersymmetric gauge theory, 
in eq. (2.19). The correspondence is easily generalized when several stacks of branes are 
present. Explicitly, in a construction with P number of stacks of branes, with number of 
branes being Ui for the i'th stack, the flux (for a given target space component {ij) ) takes 
a form: 



F- = 



F'^I 
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(A.7) 



with /„. being the rii-dimensional identity matrix and we have hidden the ij indices in 
the RHS of eq. (A.7) in constants that are all integrally quantized, as given earher 
exphcitly in eqs. (2.13) and (2.14). The corresponding gauge potentials will also then 
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have a block diagonal structure: 

(a]i 



(A.8) 



Now, in order to understand the wavefunctions associated with chiral fermion bifun- 
damentals, in such a representation of the brane stacks, we consider the flux matrix Fjj 
in eq. (A. 7) and gauge potential in eq. (A.8) with only two blocks {P — 2). The chiral 
fermion bilinears between stack-a and stack-6 are then represented by: 



ab 



Xab 



(A.9) 



with Cn^, Cn^ being constant matrices of dimensions Ua and respectively. We can easily 
derive the equation satisfied by the various Dirac components, as given in eq. (A. 5), for 
Xab such that ipab Satisfies the Dirac equation (A. 6). We obtain: 

drxl + 92X1 + - + (A' - A'),,xl = 0, 



d2xl - dixl + (A' - A'),,xl - {A' - A'),al = 0, 
dix'- + d2x'- + (A' - A').,x'- + {A' - A'),,x'- = 0, 
d2X- - diX- + {A' - A'),,x'- - {A' - A'),,x'- = 0, 



(A. 10) 



with subscript a, h being dropped from Xab to make the expressions simpler. We will, 
however, restore the indices at a later stage while evaluating the overlap of three such 
wave functions from different intersections. In particular, for the chiral components, xX 
equations reduce to: 



d,x\ + {A' - A^),,x\ = 0, 
d2x\ + {A^ - A^),,x\ = 0. 



(A.ll) 



The generalization of eq. (A.ll) to the case is straightforward and can be written 



as: 



Ax? = hxt + (^' - A^)^-iXt = 0, (i = 1, 2, 3). 
Eq. (A. 12) matches with eq. (4.65) of [5] for Vt = il^, with the identification: 



(A^-A%^|([N.(z + C)].(Iml2)-i)^, 



(A.12) 



(A.13) 
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with C being the complex constants representing the Wilson lines and is the difference 
of fluxes between the two stacks a and h (see eq. (2.15)), having constant fluxes and 
F^, giving the fermion bilinears in the representation (^1,72,2). 

Such a solution for eq. (A. 12) and (A. 13) is given in [5] for arbitrary complex structure 
Q, by the basis elements: 



with general deflnition of Riemann theta function: 



3 




fN.£N.ri) 



(A. 14) 



a 
b 



_ "y ^ g7r(m+a).n.(m+o)g27ri(m+a).(i7+6) 



(A.15) 



and N satisfying the constraints given in eqs. (2.17) as well as: 



(A.I6) 



implying that j.N is an n-dimensional vector with integer entries. Also, the normalization 
factor J\f in eq. (A. 14) is given by: 



Af= (2"|detN|.det(/mf7))4 {Vol{T^'')) 
Then wavef unctions satisfy the orthonormality relations: 



I 

Jrp2n 



(A. 17) 



(A.18) 



These results arc useful in determining the interaction terms in Section 4. However, to 
have well-defined wavef unctions, N's must satisfy the Riemann conditions given in eq. 
(2.17). 



B More information on fluxes 

In general, the (1, 1) form flux F^izJ given by a hermitian matrix in eq. (2.6) is constrained 
by two equations (2.7) and (2.8) which mix the matrix components p^^^ Pyy and p^y for 
general Q. However, for a canonical complex structure, corresponding to orthogonal tori, 
the constraints simplify and are written in the matrix form: 

Pxx=Pyy, Ply=Pxy (B.l) 
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Fluxes of such types have been used in [27] for constructing an SU (5) GUT with stabihzed 
moduh and in Section 6 we apply the Yukawa couplings computation results to show the 
mass generation for extra non-chiral states in the model of [27]. In this case, the (1, 1) 
form flux Fzizi-i for (Jl = ^73), reduces to: 

Fzizi = ^(Pxy - iPxx) (B.2) 
Explicitly, the hermitian flux matrix F in eq. (A. 7) is given as: 

^ Px^y^ Px^y'^ + iPx^x'^ Px^yS + iPx^x^^ 

F = Pa;i3/2 - iPx^x"^ Px^y^ Px^y^ + iPx^x^ ■ (B-S) 

\Px'^y^ ^Px^x^ Px^y^ ^Px^x^ Px'-^y'-^ / 

For magnetized branes in [12,27], we used the quantization rule for p's: 

_ TTT'x^y^ _ f^x^xi _ f^x^yi /-n a\ 

Px'-y^ „j ~f-i Px'-x^ — „i ~f-i Px'-y^ ~ri IJ; [^•^) 

where rrixiyj, rrixixj, iriyiyj are the first Chern numbers along the corresponding 2-cycles 
and n^' etc. are the wrapping numbers along the 1-cycles x\ y\ However, for the 
model [27], we have used only integral fluxes corresponding to n^' = n^' = 1. 

An additional modification comes when nonzero NS-NS 5-field background is turned 
on along some 2-cycle. In this case, the first Chern number along the particular 2-cycle 
(for n^' — n^' = 1) is shifted by: 

In the model that we discussed in [27], we turn on nonzero NS-NS S-field, {B — |), 
along the 2-cycles diagonally in the three T^'s. Resulting fluxes are then half-integral. 
However, as already mentioned earlier, in writing the wavefunctions of chiral fermions Xab 
in bifundamentals, the relevant quantities arc the difference of fluxes in the two stacks, 
or the two diagonal blocks in the gauge theory picture. In addition to the Z^-branes, an 
orientifold model also contains image D-branes with fluxes of opposite signature than the 
ones present in the original brane. In such cases, the corresponding wavefunctions Xab* 
will obey similar equations as that of Xab, but with the addition of the gauge potentials 
A" -|- rather than their difference as in eq. (A. 12). The relevant matrix N which will 
now be the addition of fluxes in the two stacks, rather than their difference, will once again 
be integral. 
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We also learnt from the second equation in (2.17) that (N./mfi) is a symmetric matrix. 
However, as explained in eqs. (2.6) in the general situation and in (B.2) for Q = i/a, fluxes 
are in general hermitian when components of all types: p^x, Pyy and p^y are present. 



C Fluxes for the stacks Ui ,A, B, Oi, . . . , O 



8 



In this Appendix, wc write all the fluxes in the complex coordinate basis {z, z) with 
z = X + iy iov our GUT model in [27] and used in Section 6 for the non-chiral mass 
generation. 



= "2 {dzi dz2 dzs^ 



(^dzi dz2 dz^j 



= -- (^dzi dz2 dzs^ 



(-1 \ 


(dzi\ 


1 

2 


dZ2 


\ y 


\dz3J 


/-^ \ 

2 


( dfA 


3 
2 


dZ2 


V y 


\dz3J 


A 4 3 \ 


[dzi\ 


4 i 1 


dZ2 


V 1 -y 


\dz3J 



(C.l) 



(C.2) 



(C.3) 



= -- 
2 



^03 ^ 



(dzi 



2 







(1 


4 


-3\ 




fdzA 


dZ2 


dZ3^ 


4 


1 
2 


-1 




dZ2 








-1 


-y 




\dz3J 






2 


-4 


-3i\ 




[dzi\ 


dZ2 


dZ3^ 


-4 


1 

2 


i 




dZ2 






i^3i 


— i 






XdzsJ 






(I 


-4 


2,1 \ 




(dzA 


dZ2 


dZ3^ 


-4 


1 

2 


—i 




dZ2 






^-3i 


i 


-y 




\dz3j 



(C.4) 



(C.5) 



(C.6) 
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(^dzi dz2 dzsj 



/_25 
2 

2i 



-2i 



\ 


(dzi\ 




dZ2 


/ 


\dz3J 



(C.7) 



pOe ^ 



dzi dz2 dz3 



/_25 

2 

2i 
—i 



-2i 

1 
2 

-1 



i \ (dzi\ 
—1 (i^2 



(C.8) 



-F*^'' = -- (^d^i dz2 dz3^ 



/-'4 2i -l\ 



-2i \ 



V-1 



(dzA 

dz2 
\dz3J 



(C.9) 



dz2 dzs 



/-f 2i l\ fdz,\ 



dZ2 

\dz3J 



-2z I 



\ 1 



/ 295 
' 2 



) 



dz2 
\dz3J 



(C.IO) 



(C.ll) 



-F^ = -- dz2 dz3^ 



/3 

2 



v 



33 
2 



17 



\dz3J 



(C.12) 



Using the above fluxes, one can find out the chiral multiplets in the model. This has 
been done for the brane intersections involving stacks - U^, Ui. A computation of the 
chiral fermion multiplicities on the intersections Oi — Oj and Oi — O* ,for i,j = 1, ... 8, 
implies the existence of following fields in the non-chiral spectrum of the model. They are: 
(0t^^ 40), (0ji_^% <l>'^f^ 84), (0^L^^ 84), (0jL«% 20), (0^L^^ <l>%^\ 

49), (0^L^^ 6), (0?L^«, 14), (0?i°% ^ 84), (0^i^^ ^ef ^ 84), (0^i°^ ^ef ^ 49), 

(0?_^«, 20), (0f_°^ 14 ), (0f_««, 6), (0?_«% % 40), (0?!.°% 14), (0j!.O6^ g), (0?_^^ 

20), (</.j30% </.ef % 49), (</.?i^% 6), (0?_^^ 14), (</.;i^^ </.ef ^ 49), (0ji^«, 20), (0?_^% 
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% 8), (f_f\ 20), », 20), (02!5'% % 20), «, 20) 

8), ^ 59), (^eL''^ 33), 33), % 86), (</.eL°«, 10), ^ 

24), (^JV''^ 52), {ct-K 33), (^ef^ 33), {ct-\ 10), 86), % 52), (0jf « 

24), 59), {ct>lf\ 52), 24), (0?/ % 86), {ct°\ 10), (^Jf ^ 24), {ct>°f^ 

52), (</.e!.°^ 10), (0^f 86), 41), % 23), % 23), % 23), « 

23), {4>-'-\ 41). (C.13) 

As a result of a similar analysis for the remaining stacks A and we have also the 
following fields: 

0?/, 149), (0?/*, 146), 51), 16), 0^-^, 

149), 0^^, 45), (0^?, 0^^, 2336), 18), (0^1.^*, 292), (0^^, 149). 

(C.14) 
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